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ABSTRACT 




Title of Dissertation: Charge States of Low Energy Ions from 

the Sun 


Richard Kevin Sciambi, Doctor of Philosophy, 1975 

Dissertation directed by: Dr. George Gloeckler 

Associate Professor, Department of 
Physics and Astronomy 


Measurements of the ionization states and energy spectra 
of carbon, oxygen, and (in one event) iron accelerated in ten 
solar flare particle events are reported, for energies between 
15 keV per nucleon and 600 keV per nucleon. These observations 
were made using the University of Maryland Electrostatic Energy- 
Charge Analyzer aboard the IMP-8 earth-orbiting spacecraft. 

1 

The ionization states are found to be remarkably constant from 
flare to flare, despite great variations in other event 
parameters (such as He/[C+N+0] abundance ratio). The mean 
ionization state for carbon Is 5-7, for oxygen 6.2, and for iron 
11-7 » values which are similar to the respective ionization states 
in the solar wind. The time profile of the He/[C+B+0] ratio is 
examined, and It is found that the ratio is small early in the 
event, and increases with time. The energy spectra of the 
medium ions show a flattening below 100 keV per nucleon, which is 
highly correlated with event size as measured by the "event averaged" 
flux of 130-220 keV protons. 

The atomic charge states are shown to be consistent with charge 


equilibrium established with an electron gas at a temperature of 



D s 

1.5 x 10 °K - the temperature of the quiet corona. In order to 
reconcile the observed ionization states with conditions currently 
believed to prevail at the particle acceleration site, it is shown 
that some form of particle trapping in regions of quiet coronal 
temperatures is required after acceleration. The effect of 
ionization energy loss on the particle spectra during this trapping 
is discussed, and it is found that a spectral "kink" is expected. 

It is demonstrated that the observed spectral flattening is 
consistent with adiabatic deceleration of such "kinked" spectra. 

The required value of the interplanetary diffusion coefficient is 
computed, and is shown to be consistent with the observations. 
Finally, it is shown that both the time profile of the He/(CN0) 
abundance ratio and its variation from event to event may be 
explained by the model of Nakada [l 969] in terms of variable source 
region enrichment due to coronal thermal gradients. 
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I. INTRODUCTION 


The first evidence of the presence of nuclei heavier than 
helium in solar cosmic rays was found "by Fichtel and Guss [1961] 
in the 3 September i960 solar particle event. The observation was 
made with a nuclear emulsion track detector, lifted to a high 
altitude by a sounding rocket as part of the SPICE (Solar Particle 
Intensity and Composition Experiment) program. The first evidence 
of iron group nuclei did not come until the 2 September 1966 event 
[Bertsch, Fichtel, and Reames, 1969]- 

The detectors used in these experiments were track detectors 
flown on rockets or balloons , and hence were hampered by short 
exposure times (typically 5 minutes), low time resolution, and high 
threshold energies. In spite of thes.e drawbacks, however, the 
experimental results provided interesting constraints on the 
acceleration and propagation process at high energies. These first 
measurements showed, for example, that the relative abundances of 
the energetic solar cosmic ray medium and heavy nuclei at the same 
nuclear charge to mass ratio, particularly helium, carbon, and 
oxygen, were approximately constant from event to event [Biswas and 
Fichtel, 1965] • In spite of the extreme variations from one event 
to another of the proton to helium ratio, which showed changes of 
over an order of magnitude among events (from 10 to 100 MeV per 
nucleon), [Sakurai, 197**] , the constancy of the abundances of heavier 
elements was not unexpected, because of the similar rigidities of 
helium and fully stripped heavies. The differential energy per 
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nucleon spectral shape of the medium-group (6 <_ Z <_ 10) nuclei 
was found to be the same as that of helium, within uncertainties, 
above 10 MeV per nucleon [Biswas and Fichtel, 196U]. 

In addition, the abundances of the nuclear charge 6-26 group 

seemed to mirror the solar photospheric abundances as determined 

by spectroscopic measurements. It is believed that these medium 

and heavy nuclei in the solar cosmic particles are a good sample of 

the relative abundances at the sun, and therefore are a more direct 

means of sampling these abundances than is otherwise available. 

In fact, Bertsch et al [1972] and Crawford et_ al_ [1975] have 

suggested revisions in the relative solar abundances of neon, argon, 

and sulfur, 6ased on their measurements of these nuclei in solar 
l 

radiation above 15 MeV per nucleon. Another important consequence 
of this work has been the determination of the relative abundance 
of helium, which cannot be measured spectroscopically [Aller, 1971]* 
The problem of medium and heavy nuclei in solar cosmic radiation 
thus seems to be well understood at energies greater than about 
10 MeV per nucleon. 

Recently, -the measurements of medium and heavy group nuclei 
in solar flares have been extended to lower energies (less than 
10 MeV'per nucleon) by several groups, using both passive detector 
systems (i.e., track detectors and nuclear emulsions) and active 
systems (telescopes, etc.) on balloons, rockets, and satellites. 

The most interesting result of these low energy measurements is the 
indication of an enhancement, over photospheric values, of the medium 
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and heavy nuclei relative to helium — an enhancement which is 
variable from flare to flare. 

Armstrong and Krimigis [l97l] summarized the He/CNO data for 
thirty-five events during the period 1967-1968, and were the 
first to show that the composition at low energies (0.5-2. 5 MeV per 
nucleon) was different from that at higher energies. Armstrong 
et_ al^ [1972] compared their satellite data for this ratio with the 
SPICE emulsion results from I96O-I969. The satellite data have an 
average He/CNO ratio of 27 ± 9 while the rocket data average 58 ± 5 
for the energy range 12 to 100 MeV per nucleon. Both energy and 
event variations seem likely. Teegarden et al_ [1973] reported 
He/0 ratios at 13.5-47 MeV/nucleon which differ significantly 
between the April 1971 and September 1971 events. For the latter 
event they obtained a He/CNO ratio of 26 ± 2. They also observe 
Fe/0 ratios for these two events which differ by a factor of six, 
but which are statistically consistent with a constant value. 
Mogro-Campero and Simpson [1972], however, derived Fe/O ratios 
above 5 MeV per nucleon for seven events and showed significant 
variability — as much as a factor of 20. Their results were 
obtained by integrating over the major portion of each event, and 
using the data from IMP-4 and IMP-5 to extrapolate their 0G0-5 
data for oxygen from l4 MeV per nucleon down to 5 MeV per nucleon. 
The overall averages from all seven events indicated that solar 
species above silicon were preferentially enhanced relative to the 
photospheric values. The enhancement factor for iron is between 


10 and 20. 



Crawford et^ al [1972] found no evidence of composition 
variations between the January and September 1971 events, using plastic 
detectors aboard SPICE rockets. They did observe, however, that 
the Fe/Si ratio decreased with energy over the range 1.5-50 MeV per 
nucleon. The overall change is by about a factor of U. 

Price et al [1971] analyzed tracks left in the Surveyor 3 TV 
camera filter, which was exposed on the moon for about 2.6 years. 

They compared the iron spectrum to a helium spectrum obtained 
from satellite data integrated over the same time period, in the 
energy range 6-10 MeV per nucleon. The He/Fe ratio inferred from 
their results ranges from about 230 at 6 MeV per nucleon to about 
560 at 10 MeV per nucleon. The latter value agrees well with that 
obtained at higher energies. 

Van Allen et al [ 197 ^ ] have recently reported results which 
indicate that the He/medium ratio increases with time over the 
course of a single flare. The implication is that the medium 
group nuclei from a flare arrive at the earth before the helium. 

Many attempts have been made at explaining the observed 
composition enhancements [Price et_ al , 1971; Cartwright and Mogro- 
Campero, 1972; Van Allen et al ,197^ ] • Most of these mechanisms involve 
incomplete stripping of the heavy ion's electrons, followed by some 
form of Fermi acceleration. Thus, when the particles reach the 
acceleration region, the partially ionized heavy nuclei are 
preferentially accelerated, since more of these have a rigidity 
greater than the cut-off rigidity for acceleration [Wentzel, 1965 ] 
at a given energy per nucleon. 
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The ionization states hre also a sensitive measure of the 
temperature at which a particle has "been accelerated or stored, the 
amount of matter it has travelled through, and its rigidity (since 

rigidity varies as mass/ ionic charge - e.g. at 100 keV per nucleon, 

+12 +26 
Fe has a rigidity 'of 6k MV, while Fe has rigidity of -30 MV and 

He +2 has 27 MV). 

In view of the importance of the role of a particle's rigidity, 
in particular, in acceleration and propagation processes, it seems 
essential to measure the ionic charge states of solar particles and 
obtain accurate values for their charge to mass ratios as a function of 
energy. Indirect measurements have been made in the past. Sullivan 
and Price [1973], using observed spectral turn-downs on a rocket 

+4 

flight, reported charge states for iron at 2 MeV per nucleon of 22_^, 

based on calculated cut-off rigidities. Braddy et al [1973] have 

presented evidence, based on spectral shapes, that iron should be 

completely stripped while being accelerated. 

The University of Maryland electrostatic energy-charge 

analyzer experiments on board the earth-orbiting satellites of 

IMP-7 and IMP-8 are the first devices designed to measure the 

ionic charge states of low energy particles , down to less than 

> 

40 keV per charge. Results from these experiments, which have been 
presented earlier, indicated that low energy carbon and oxygen 
are nearly fully stripped [Gloeckler et_ al, 1973; Gloeckler et_ al , 

197 1 *; Sciambi et_ al , 197^*, Gloeckler et al, 1975a]- This work 
represents an attempt at a systematic, statistical study of the 
charge states, energy spectra, and abundances of low energy 
.carbon and oxygen from the sun. 



Chapter 2- of this work, supplemented by Appendix A, gives 
a complete description of the IMP-8 experiment . 

In Chapter 3, data from ten solar particle events observed 
by the IMP-8 experiment are presented. These data contain several 
new features . The ionic charge states of carbon and oxygen are 
found to be remarkably consistent from event to event, in spite 
of wide variations in other event characteristics such as alpha/CKO 
ratios, event size, etc. Spectral index data from the lowest 
energy CNO channels indicate a spectral "turn-down" not attributable 
to velocity dispersion which is highly correlated with the low 
energy proton flux. A time variation of the alpha/ CKO ratio 
throughout an event is found which is consistent with the 
observations of. Van Allen et_ al [197*+ ] • Measurements of the 
ionic charge states of iron at less than 10 keV per nucleon are 
made for the "iron-rich" event reported by Gloeckler et^ & 1 _ [1975b] , 
which indicate that the mean ionic charge state for low energy iron 
is about 10 to 12. 

In Chapter U the charge state results for carbon, oxygen, and 
iron are shown to be consistent with ionization state equilibrium 
established in the quiet corona, at an electron temperature of 1.5 
x 10 6 °K. If the particles are accelerated in regions of high 
coronal temperatures (T > 10 °K) [Cartwright and Mogro-Campero , 

1972, McGuire et_ al_, 1975] or stripped during acceleration, then 
some form of post-acceleration trapping is required. It is demo- 
strated that this trapping, combined with adiabatic deceleration, can 
account for the features observed in the heavy ion spectra. 
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II. INSTRUMENTATION 

The University of Maryland experiment on board the IMP-8 
spacecraft is fully described in Appendix A. The characteristics 
are briefly repeated here for the convenience of the reader. The 
experiment is an electrostatic deflection sensor. An ion entering 
the instrument has its entrance trajectory defined by a collimator 
(see figure A-l). The particle then enters one of three electro- 
static field regions where its path is deflected by an amount 
inversely proportional to its kinetic energy per charge. On exit 
from the deflection cavity, the ion strikes one of an array of 
rectangular solid state detectors, each of which covers a well- 
defined deflection, and hence energy per charge, range. Pulse-height 
analysis of the detector signal allows a direct measurement of 
the particle's total kinetic energy, from which its charge may be 
determined. In addition, by using multiple energy thresholds on 
signals of several detectors, counting rates having well-defined 
energy limits and hence charge responses are generated with ti m e 
resolutions as short as five seconds. Several of these rates are 

sectored into four directional quandrants in the ecliptic plane, 

> 

providing data on particle anisotropies. A plastic scintillator 
anticoincidence cup and an array of solid-state anticoincidence 
detectors placed behind the detector array prevent analysis of 
penetrating particles. 

Background correction of a detector's pulse-height histogram 
is performed by subtracting a normalized histogram obtained from 
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the same detector during quiet times, when the counting rates are 
at their lowest values. During these quiet times, analysis of 
pulse-height histograms indicates that the detectors respond 
primarily to the secondary electron background produced in the 
spacecraft . 

The relative abundance of individual charge states may be 
determined by fitting computed instrument response curves to the 
observed pulse-height distributions, as described in Appendix C. 

A full description of the instrument response characteristics is 
given in Appendix B. The sensor is capable of fully resolving 
ionic charge states 1 to 4 for elements up to oxygen. In addition, 
charge states 5 to 8 are separated sufficiently so that, by means 
of a simple iteration procedure, abundances of these individual 
states may also be determined. (For example, the charge-state 
resolution of detector P3 is plotted in figure 2-1.) Complete 
charge state identification for carbon and oxygen down to about 
40 keV per charge is thus obtained. Of course, the individual 
elements may not be resolved in this manner since only the ionic 
charge (as opposed to nuclear charge) is determined. In many cases, 
however, the ions may be completely identified (i.e. , ionic charge 
and nuclear charge) by using the simultaneous nuclear abundance 
information provided by the Ultra Low Energy Telescope (ULET) 
[Hovestadt and Vollmer, 1971] of the Maryland/Max-Planek-Institut 
experiment at somewhat higher energies. 
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ULET is dE/dx vs total energy telescope with a geometrical 
2 ' 

factor of .53 cm sterad. The sensor is pictured in figure 2-2. 

The front, AE element (Dl) is a thin-window flow-through proportional 

counter, with a total mass of a, 300 pgm/cm . The rear, E element 

(D2) is a silicon surface barrier detector. A plastic scintillator 

anticoincidence cup (A) surrounds the detector assembly. The low 

mass of the AE element allows two parameter analysis of incoming 

particles down to < 100 keV per nucleon for iron. Different 

elemental tracks can be reliably (i.e., automatically) separated 

down to about 600 keV per nucleon. Every three hours ULET changes 

its analysis threshold, so that it alternately analyzes all incoming 

particles, and only those particles heavier than helium. ULET 

shares the electrostatic energy-charge analyzer's pulse height 

3 1 

analyzers and telemetry on a basis, and both instruments 

are contained in the same spacecraft module. 

IMP -8 (Explorer 50) was launched from Cape Canaveral , Florida 
on 1973 October 26 at 0226 GMT. The spacecraft achieved an 
eiliptical orbit (shown infigure 2 - 3 ) with a mean apogee of about 
^5 earth radii and a perigee of a, 28 RE, inclined 29° to the 
ecliptic. The University of Maryland experiment was turned on 
on October 30 and has operated normally since that time. 
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III. OBSERVATIONS 

The majority of the ten events used for the study were of small 
or moderate size. The reason for this selection is that during periods 
of large high-energy proton fluxes, the shape of the background histogram 
for a detector changes, making correction of the histogram for background 
difficult. In addition, in an attempt to exclude from the analysis 
particles of magneto spheric origin [Fan et al, 1975] , only events 
observed when the spacecraft was outside the earth's bowshock- 
magnetotail region (SEV angles 2k0° - 120° in Figure 2-3) were used. 
Characteristics of the events analyzed are shown in Table 3-1. 

The PHA Data 

For each event, pulse-height histograms were accumulated over 
the major portion of the active period'. Each histogram was corrected 
for background contributions by subtracting a suitably normalized 
quiet time histogram from the same detector (see Appendix C). 

Figure 3-la shows a typical uhcorrected puLse-height spectrum 
(histogram) obtained during the solar flare particle event of 197 ^ 

June 22. The measured, smoothed background distribution is shown 
as a dashed curve, and the histogram after background correction 
is shown in Figure 3-lb. This histogram- was obtained from IMP-8 
detector P3, which responds to ions in the energy range 130 to 220 
keV per charge. The triple-peaked shape of the histogram is 
characteristic of every event studied to date. The two left- han d 
peaks are produced by protons and alpha particles respectively, 
while the right-hand peak is due to particles having ionic charges 
5-8 for carbon and oxygen, and 8 _< Q 12 for iron. The large dip 
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between the alpha peak and the Q=5-8 peak is statistically consistent 
with a complete absence of ionic charges 3 and 4. In none of the 
events examined has there ever been any evidence for the presence of 
these two charge states . At this point, computed instrument response 
curves for ions of various charges were fit to the data in order 
to determine ionic abundances. The first stage of the fitting of ' 
charge states to the PHA data was to determine the abundances of 
ions of charge 5, 6, 7, and 8, regardless of nuclear charge. In 
order to do this, all ions with charge states 5 and 6 were assumed 
to he carbon and all ions with charge states 7 and 8 were assumed 
to be oxygen. (Some assumption of the nuclear charge is necessary 
to correct for the nuclear defect [Lin and Campbell, 1973; see 
Appendix B].) Next, to determine the amount of oxygen present with 
an ionic charge of 6, elemental abundance data for the same time 
period were obtained from ULET at _> 600 keV per nucleon for the ratio 
of carbon to oxygen. This ratio was assumed to remain constant 
down to energies appropriate for the deflection sensor, and was used 
as an input to the program MSPCT (see Appendix C) in order to determine 
the relative abundances of C and 0 . After analysis of the iron 
rich 197^ May l4 event, it was found necessary to correct the charge 
state data from most other events for iron contamination. 

The Ionic Charge States of Iron 

The 1974 May lU event is part of the "iron rich" event reported 
by Gloeckler et al [ 1975b]. They find during this time period a 
sustained emission of particles with a composition richer in iron 
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than in oxygen. The small event on May 13 produced sufficient fluxes 

for charge state analyses to he performed. During the time period 

for which PHA histograms were accumulated, ULET data at 600 keV per 

nucleon indicate that the iron/oxygen ratio is about 2.5, and that the 

medium-group spectral index is about 2. If the composition remains 

so enriched down to the energies seen by the deflection sensor, 

then most of the heavier ions seen by the electrostatic sensor 

should be iron. The solid histogram in Figures 3-2a and 3-2b shows 

the observed pulse-height distribution from detector P3 during this 

period. Note the presence of a large peak due to heavy ions centered 

at about 1 MeV. By using the ULET abundance data for C/0 and Fe/O, 

it is possible to generate predicted response histograms based on various 

ionic charges for -iron. For example, if iron were fully stripped, i.e., 

+26 

all iron was Fe , the observed histogram should have taken the form of 
the dashed curve in Figure 3-2b. Clearly the required Fe + ^ peak (the 
rightmost dashed peak) is not present in the observations, and the 
measured oxygen flux is above the predicted value. Thus it must be con- 
cluded that iron is only partially stripped at 150 keV per charge; 
in order for Fe ions to appear in the 1 MeV peak, the mean ionic charge 
must be about 10. In the dashed curve of Figure 3-2a, iron charge states 
of 8 to 12 have been fit to the data. The measured charge state distri- 
bution then is similar to that of the solar wind [Holzer and Axford, 1970] ; 
our measured relative abundances of the +10, +11, +12, and +13 ionization 
states for iron are 15%, 35%, 30%, and 20% respectively. The fit to the 
data in this case is extremely good. Using these measured charge states, 
it is found that the agreement between the deflection sensor and ULET for 
the iron fluxes is excellent (Figure 3-3). So by requiring that the spectra 
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from the electrostatic sensor and the ULET join smoothly, the 
first direct measurement of the charge states of low energy iron 
has been made [Gloeekler et al , 1975e]. 

Iron Contamination of the Carbon and Oxygen Peaks 

The fact that iron is present in these low charge states (i.e. 
overlapping the carbon and oxygen peaks) makes it necessary to 
correct the charge states of carbon and oxygen determined in other 
events, even though the measured iron abundances are small. The 
reason is that for a steep spectrum even a small abundance of 
partially stripped iron can make a large contribution to the 
carbon-oxygen peak. 

The deflection sensor measures the kinetic energy per charge 
of an incident ion. Abundance ratios are usually calculated at 
equal energy per nucleon . Suppose that the incident ions are species 
■1 with charge Q x at kinetic energy Eq, and species 2 with charge 
Q2 a "t kinetic energy E^. Then in a given detector,- 


E, 


E„ 


_1 _ 

Qq ^ 2 


Defining the energy per nucleon £ - E /A 

n n' n 


3 , 


e 2 A g 

q 2 


or. 




If the incident flux of particles is a power law in energy per nucleon, 

F - 

then the ratio of fluxes in the detector will be 

fi. _ ®i /fi\ Y 

F 2 a 2 \ 2 / 

= /^\" Y 
a 2 \ Q 2 A i/ 


= a. 
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where a ^ is the abundance ratio of species 1 to species 2 at equal 


energy per nucleon and n = 


_/V 2 y Y . 


\ Q 2 A 1/ 


is a spectrum dependent abundance 

+12 . 


enhancement factor. For a spectral index of y = 3, Fe is apparently 

+6 

enhanced over its actual abundance relative to 0 by a factor of 5.3. 

Thus, in all events with a steep spectrum { y greater than l) or a large 

iron abundance, the computed carbon and oxygen histograms must be 

corrected by subtracting iron in appropriate charge states. This will 

+8 

tend to reduce the abundance of 0 , as the observed iron charge states 

fall at the high end of the carbon-oxygen range. 

The Rate Data 


The pulse-height analysis data are approximately an order of magni- 
tude more sensitive to variations in cosmic ray intensity than the rate 
data. In addition, these data represent the only way to determine 
the ionic charge states of the incident particles. However, the pulse- 
height analysis is restricted by spacecraft telemetry to two 



analyses every 1.28 seconds. This low rate of accumulation means 
that the use of PHA techniques is statistically not feasible for ' 
short time periods (i.e., less than 6-12 hours when the EECA has 
3 A of the telemetry time). Because of this limitation, the PHA 
data are insensitive to short term variations of event characteristics; 
hence the results based on PHA analysis represent an average over 
the event. It becomes necessary therefore to use the rate data in 
addition to the PHA data in order to verify event characteristics, 
to obtain time profiles of the helium/medium ratio, and to examine 
velocity dispersion, anisotropies, etc. 

The counting rate data are corrected for background very simply 
using the P5 detector rate, which counts only background. Since the 
background (due to secondaries and penetrating particles which are not 
ant icoinci dene ed) should be more-or-less omnidirectional, the ratio 
of the background counting rates in two detectors should be a function 
of their thresholds and the ratio of their uncollimated geometrical 
factors, or their active areas. Thus, a given rate R^ is corrected for 
background by subtracting a well-determined fraction of the P5 rate R^_. 

In addition to plotting particle fluxes determined from the 
rate data, certain ratios of fluxes are also plotted for the events. 

For example, the ratio of the P2 flux to the P2" flux gives the 
helium/medium-group ratio. The spectral index between any two rate 
channels can be determined from 

Y = -[log(J 2 /J 1 )]/[log(E 2 /E 1 )3 

where and E^ are the measured background-corrected fluxes and 
center point energies of the rate channels. 
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So, by using both the rate data and the pulse height data, a 
reasonably complete picture of the event emerges. 

The Events 

The first event analyzed was the 1974 July 4 increase. This 
particle event was the result of a very complicated series of flares 
involving three sudden commencements (abrupt increases in the earth's 
magnetic field due to the arrival at the bowshock of a high-speed 
plasma stream) at least one of which was associated with an inter- 
planetary shock. During this event the highest fluxes of any event 
since the launch of the IMP-8 spacecraft were observed. Rate data for 
15-600 keV per nucleon ions heavier than helium indicate a strong 
bend-down in the energy spectra below 100 keV per nucleon, even late in 
the event, which cannot be attributed to velocity dispersion. These rate 
data are shown for the entire event in Figure 3-4, where the fluxes 
at 40, 80, and 150 keV per charge, corrected for background, are plotted. 
The two lowest energy fluxes track each other throughout the event 
at the same level. Rate data for protons and helium show no such 
bending within the energy range of the detectors, (it should be noted 
however that the low energy limits for these two species are consider- 
ably higher than for CNO.) The helium/CNO ratio at, 80 keV per charge 
is plotted in Figure 3-5* This ratio starts low at the beginning of 
the first particle increase on day 185 . The ratio grows with time, 
reaching its local maximum late in the first rate peak. It then drops 
again as the second large increase begins, and once more grows larger 
with time during the particle increase. 
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Pulse-height analysis data support these observations, and indicate 
that carbon and oxygen are present only in ionic charges >5. The 
event averaged energy spectra for helium and CNO are shown in Figures 3-6 
and 3-7 respectively. After correcting for iron contamination, we 
find that at 150 keV per charge the abundances of charge states. 5 
and 6 for carbon are 28$ and 72$ respectively, and that the abundances 
of oxygen charge states 6, 7, and 8 are 68$, 27$, and 5$ respectively 
at the same energy. The abundances for these ions at 500 keV per charge 
are 25$, 75$, 70$, 24$, and 6$ respectively: the differences are not 

statistically significant. The mean charge states of carbon and oxygen 
are plotted as a function of energy in Figure 3-8 — we see little variation. 

The PHA data were also divided into two parts, one early in the 
event, and the other late (the dividing time was 0900 on day 187, July 6). 
Comparison of ionization states derived from these two parts showed no 
significant differences: the ionic charge states appear to be constant 

throughout the event. 

It is instructive to examine a simpler event — one in which local 

effects are not likely to be as important as in the July 4 flare. In 

Figure 3-9, for example, are shown proton, helium, and CKO spectra 

ll ll 

averaged over the time period 1974 June 10 0 to 1974 June 12 0 . 

This event, associated with a small flare on June 8, is typical in 
that the CNO spectrum flattens at lower energies. An unusual feature 
of this event (as compared with the others examined) is that the 
helium also seems to show a flattening; in fact, this is the only 
case in which such an effect has been seen in the alphas in the energy 


range of the instrument. 
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Charge state data from pulse-height analysis show that the 
relative abundances C + ^/C, C + ^/C, 0 + ^/0, 0 + ^/0, and 0 + ^/0 are 
n%i 83$, 82$, 1 6 $, and <2% respectively at 150 keV per charge. 

Again, at high energies (500 keV per charge) there are no signi- 
ficant differences in these fractions. Mean charge states for carbon 
and oxygen are plotted vs energy in Figure 3-10. 

Figure 3-11 shows the time variation of the spectral index y 
(J ^ E Y ) for heavy ions between 20 and 40 keV per nucleon during the 
June 10 event. The effect of velocity dispersion is evident up to 
the beginning of day l6l, at which time the spectral index settles 
down to a low, more-or-less constant value and remains at that value 
to the end of the event. The bend-over in the spectrum (low spectral 
index y) cannot, (therefore, be attributed to simple velocity dispersion. 

Using the arrival times of the low energy heavy ions of various 

energies, it is possible to calculate their diffusion coefficients, 

k h , based on travel distances (along the average spiral magnetic field) 

of 1.2 AU for this event. The increase in FI is seen about 37 hours 

after the flare, while that in P2" is seen about 4.5 hours earlier. 

The P 6 increase is seen about 25 hours after the flare. Inserting 

these times into the relation for diffusive propagation 

R 2 

K ll = £7 

20 20 2 

values of k jj of between 4 x 10 and 6 x 10 cm per second are obtained 
for ions with energies between 20 and 250 keV per nucleon. The very 
small {<3%) residual anisotropies which are observed after correction 
for the Compton-Gettirig effect [ipavich, 1974] are consistent with 
this low value of the diffusion coefficient. These values are similar to 
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those obtained for other events at low energies by Lin et_ al [1974] 
and Lupton and Stone [1973]. 

In Figure 3-12 is plotted the helium/medium ratio as a function 
of time for this event. Again, the ratio tends to increase as the 
event progresses, as though the heavier ions arrive before the alphas. 

Figures 3-13 and 3-l4 show proton, helium, and CKO spectra 
from two other events in the study, 1974 June 22 and 1974 Kovember 8. 

The June event shows no turndown in the spectra, even at the lowest 
energies. The Kovember event shows a large bendover in the Q=5-8 
group spectrum, similar to that observed in the July 4 and June 10 
events. However, rate data indicates that there may be some local 
effects present in this event. In both these events, as well as in 
all the others examined, the mean charge state of carbon is about 5-7, 
and of oxygen is 6.2. This constancy of the ionic charge states 
stands out, despite great variations in event parameters^ 

Ionization state composition data for all ten events are shown 

| r | ^ 

in Table 3-II. For each event, the relative abundances C /C, C /C, 

0 /0, 0 /0, and 0 / 0 have been averaged over energy, as there were 

no statistically significant variations of the charge states with 
energy. The abundance of each of these ions averaged over all events 
is shown in the last line. We note once more the remarkable con- 
sistency of the ionic charge states from event to event. Also included 
in Table 3-II are measured ion abundances in the solar wind, and the 
expected charge state abundances in the solar corona. 

Table 3-III summarizes the spectrum data for the ten particle events, 
al 1 of which were associated with solar flares on the western disc of the 
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sun. Four of these events show an "unusually low" spectral index 
(y < l) at low energies for CKO, while, with the exception of the June 
10 event, none show a clear alpha turndown. We must caution, however, 
that since for alpha particles the low energy limit of the observations 
is 40 keV per nucleon, a factor of two higher than the lowest heavy 
ion energy, the possibility of an alpha turndown below 40 keV per 
nucleon cannot be ruled out. On the other hand, one can definitely 
say that if a turndown is present in the alpha spectra, it is not nearly 
as strong as for the case of the heavy ions, or occurs at lower energies. 

The alpha/CNO ratio at 40 keV per nucleon shows a great variability 
from event to event, but there is no correlation of CKO abundance with 
the degree of turndown in the spectra; nor is there a correlation with 
helium flux, CKO flux, or location of the flare. There is, however, 
a strong correlation of CKO spectral index Yqjjq with the low energy 
(<220 keV) proton flux averaged over the event. This interesting 
correlation is evident in Figure 3-15, where the CKO spectral index 
from U0-80 keV per charge is plotted va the "event-averaged" 130-220 
keV proton flux. The trend is clear: the larger the flare as measured 

by the low energy proton flux, the stronger the bend in the spectrum. 
Summary of Observations 

The major observations may now be summarized: l) The mean ionic 
charge states of carbon and oxygen are 5-7 and 6.2 respectively, showing 
no variation with energy from 37 keV per charge to greater than 1 MeV 
per charge, and only statistical fluctuations from event to event; 

2) The mean ionic charge for iron over the same energy range is found 
to be 11.7*0.6 (for one event); 3) In the larger events, a turndown is 
observed in the medium-group differential energy spectra below about 
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100 keV per nucleon, which seems to he well correlated with the event 
size as measured hy the "event-averaged" 130-220 keV proton flux; 
and 4) In two events a time variation in the alpha/medium group ratio 
at 80 keV per charge is found which indicates that the medium group 
nuclei arrive before the helium. In addition we see a large variability 
from event to event in the heavy ion abundances relative to helium. 
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XV. DISCUSSION 

As was pointed out in the previous section, the ionic charge 
states of the medium group nuclei measured in the events of this 
study show remarkable consistency from event to event. In this 
section the significance of the observed charge states and the 
features of the heavy ion spectra is examined, and it is shown that 
these measurements are a sensitive probe of the conditions of the 
acceleration region. 

The Charge States and Characteristic Temperature of the Ions 

Can the charge states be explained in terms of flare region 
conditions? According to current theory on x-ray emissions, during 
the flash phase of a flare, large numbers of keV electrons must be 
generated [Takakura, 1967; de Feiter, 197**; Kane, 1973, 197**]- It 
has been suggested [Braddy et al , 1973; Sakurai, 1975] that these 
electrons could strip heavy ions by impact ionization. However, 
the observed charge states, while high, do not indicate complete 
stripping, as would be expected if this mechanism were the dominant 
one. 

It may be noted from the work of Jordan [ 1969 ] and others 
[Allen and Dupree, 1967; Cox and Tucker, 1969; Lange and Scherb, 
1970] that the ionization states observed are consistent with those 
expected for ions in thermal equilibrium with an electron gas at 
a temperature of 1-2 x 10^ °K. In addition, the charge states are 
very similar to those observed in the solar wind [Bame et_ al, 1968 , 
1970], although the elemental abundances are different. It thus 
seems likely that the low energy heavy ions seen near the earth have 
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achieved charge state equilibrium in regions of the solar corona, 

from -which the solar wind orginates. In figure 4-1 are plotted, 

after Nakada [1969 ], the expected charge states vs_. temperature for 

several, elements, with our data superposed. The data, corrected 

for iron abundances in each event as measured by ULET at 600-1000 

keV per nucleon, show that both carbon and oxygen charge states fit 

well to a temperature of about 1.6 x 10 °K. How, using measured 

+4 +8 

upper limits on the relative abundances of C and 0 , both lower 

and upper limits can be set on the temperature of the region in which 

the ions achieved charge state equilibrium. Specifically, since the 

ratio C +4 /(C +4 +C +5 +C +6 ) is observed to be less than 10$, we can set 

T . > 1.0 x 10^. And since the ratio 0 + ^/( 0 + ^+0 + ^+0 + ^ ) is less 

min 

tr 

than 10$, T < 2.5 x 10 . 
max 

How sensitive are the charge states to the source region 
temperature? Kafatos and Tucker [1972] have done a. study of the 
response of the ionization equilibrium states to changing plasma 

temperatures. They find that if the region containing the particles 

7 9—3 

is Instantaneously heated to 10 °K at an electron density of 10 cm 

(typical parameters for the flash phase of a flare), the characteristic 

ionization time t^ for the heavier ions is very short. For carbon, 

for example, t^ is of the order of 10 seconds - much shorter than the 

duration of the flash phase. (By contrast, the ionization time 

g 

for carbon to be heated from lower temperatures to 10 °K at the 
same electron density is about 800 seconds.) Thus immediately after 
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heating the charge states of the medium group ions should reflect 
the higher temperatures — carbon and oxygen should be fully stripped. 

Even if the particles immediately return to a region of temperature 
and density characteristic of the quiet lower corona, there is 

Q 

considerable hysteresis in their ionizabion states. For n =10 

e 

3 6 

cm , T = 1.5 x 10 °K, the recombinat Lon time t^ for oxygen to 

return to its equilibrium charge states is of the order of 5000 
seconds. This time may be considerably lengthened, for if the 
particle moves radially outward the electron density it sees rapidly 
decreases, and the recombination time varies inversely as electron 
density. 

We thus have the problem of explaining how the ions we observe, 

which were clearly flare accelerated, can have ionization states 

which are typical of the quiet corona. The first possibility is 

that the temperature oy the acceleration region does not exceed 

6 / 

1.5 x 10 °K - the /fuiet corona temperature. 

Although there is no direct experimental evidence allowing 
determination of the site of flare acceleration of heavy particles, 
there is reason to believe from x-ray and radio burst observations 
of the flash phase that the acceleration of electrons, at least, occurs 
in the lower corona, at temperatures of 1-3 x 10 °K and densities 
greater than 10^ cm ^ [Kane, 1973; Kahler, 1973]. This is consistent 
with current flare models [Sturrock, 1973], and with observations 
of flare site temperatures [Strauss and Fapagiannis, 1971]- Thus 
the only available evidence contradicts acceleration at low temperatures. 
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The second possibility is that the particles, after acceleration 

£ 

at high temperatures, are trapped in regions of T ^ 1,5 x 10 °K 

11 -3 

for a time density product of the order of 5 x 10 sec cm . This 
trapping region may be the acceleration site itself, cooling back 
to pre-flare temperatures; alternatively the particles may escape 
the acceleration site only to be trapped elsewhere in the corona. 

Even if the first possibility is accepted and the ions are 
accelerated at low temperatures , if they have been stripped by 
flare-associated keV electrons some storage may be required to 
reestablish the characteristic charge states. 

There is evidence for trapping of flare-accelerated particles. 
Pulsations in microwave bursts may indicate temporary containment 
of energetic particles in closed mangetic loops [Newkirk, 1973]- 
Stationary type IV bursts have been interpreted as synchrotron 
emission of relativistic particles trapped in stable magnetic arches, 
precipitating into the corona — these may last for horn's or days 
[Wild and Smerd, 1972]. Stable closed loops suitable for trapping 
have been observed out to 2.5 Rq [Boischot, 197^]. The possibility 
of trapping or storage of low energy protons has, been suggested 
by Fan et al [1968], Anderson [1969], McDonald and Desai [1971] s 
Simnett and Holt [l97l] s Simnett [1971 S 1972, 1973], and Ma Sung 
et_ al_ [19753 among others. 

Spectral Shapes, Particle Storage, and Energy Loss 

If the ions are stored somewhere in the corona after acceleration. 


it is expected that their energy spectra will suffer some degradation 
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due to ionization energy losses. The presence of such a degradation 
might give evidence as to whether storage of the particles actually 
takes place. Is the observed spectral turndown produced by the 
effects of ionization energy loss during storage, or is it 
caused by propagation effects in interplanetary space? 

Verzariu and Krimigis [1972] have pointed out that by allowing 
the diffusion coefficient k to depend on particle velocity, a bend 
in the spectrum can be produced by adiabatic deceleration. But, 
since the energy change mechanism would depend only on particle 
velocity, different species observed at the same velocity should 
have experienced identical propagation effects, and thus should 
"map back" to the same source velocity. However, in the present 
case, the absence of a strong alpha turndown at the same kinetic 
energy per nucleon (or velocity) as for carbon and oxygen makes it 
unlikely that the bending is caused by a velocity-dependent k. If 
the initial spectra were the same, then the same mechanism which 
produces a CNO turndown should produce the same degree of turndown 
in the alphas. Using an identical argument, the effect of a 
rigidity-dependent diffusion coefficient can be shown to be 
unlikely, since by virtue of their high state of ionization the 
carbon and oxygen have nearly the same rigidity (within 20-30% ) as 
the alphas at the same energy per nucleon. Adiabatic deceleration 
can certainly change the shape of the spectra: The point here is 

that this deceleration should produce identical changes in the 
spectra of all species observed. 
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Having ruled out propagation effects, we now examine 

coronal trapping. Previous attempts to calculate the effect of 

storage on the spectra of flare particles have led to the conclusion 

that such storage could not explain the shape of the proton 

spectra observed at the lowest energies [Ahluwalia, 1971; Krimigis 

and Verzariu, 1971; Verzariu and Krimigis, 1972; Krimigis, 1973; 

Newkirk, 1973]. We point out, however, that in each of these cases, 

the authors have used a form of the energy loss equation which is 

invalid at low particle energies. 

In calculating the effect of ionization energy loss on the 

spectra of the particles during propagation from the acceleration 

region and/or storage in the corona, it is important to use the 

correct form of the energy loss equation for low particle energies. 

The form valid at high energies, dE/dx ^ l/E, breaks down as the 

particle velocity approaches the electron thermal velocity (v. . ). 

the 

g 

For an electron temperature of ^10 °K, this corresponds to a 

particle energy of a, 225 keV per nucleon. 

The problem of the slowing down of an ion in a fully ionized 
plasma has been treated by several authors. We use here the 

3 

results of Itikawa and Aono [1965] , who find a form of the energy 
loss equation which is valid for all values of the ion velocity. 

In the case where losses due to plasma ions are ignored, in the 
limit M_. » m^ this equation becomes: 
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u 


(4-1) 


with u = 79*5 /E./T , where E. is the ion energy in keV per nucleon, 
1 e 1 

T e is the electron temperature of the plasma in °K, is the ion 
mass in amu, x is the path length in cm, n^ is the electron density 
of the plasma (per cm' ) and erf is the normal error function: 
erf (u) = — 

the specific energy loss dE^/dx is proportional to v^ Qn or /eT, 

while in the high energy limit (v. » v,, ), dE/dx ^ l/E.. The 

^ 0,7 xon the x 

last bracketed term in (4-l) is plotted in figure 4-2. This 

function has the form predicted by Ginzberg and Syrovatskii [1964] 

for ion velocities near the plasma electron velocity. 

It is instructive to first examine the behavior of a spectrum 

of particles in two limiting cases: One in which the particle 

population is propagated directly outward from the acceleration 

2 

site to earth, passing through ^80 pgm/cm of matter, and the 
other in which the entire population is trapped for some time 
(independent of particle energy) before being released to propagate 
to earth. In the latter model the amount of matter passed through 
is dependent on particle velocity. 

Consider first the case in which all particles are propagated 
through the same amount of matter. The low energy form of the 
ionization loss equation (4-1) predicts 


/ exp(-x^) dx. In the low energy 




xon 


dE 

dx 


a/E . 
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After passage through a slab of matter, a particle of initial 
energy E^ has its energy degraded to an energy E 

E 0 = (ax + Ji) 2 ; • 


Now consider a differential spectrum of such particles. The initial 
spectrum J^(Eq) is degraded, and the new spectrum J (by conservation 


of particle number) obeys 


J 0 (E 0 )dE Q = J(E)dE 
J(E) = J„(Ej 


dE 


0 


O' O' dE 


or 


J(E) = J [(ax + /S) 2 ] ^ . 

0 Ji 


If Jq(Eq) is a power law in E^ 


- then , 


W = A V Y 


J(E) = A (ax + /e)’ 27 (ax + /e) E " 1 ^ 2 
= A (ax + /e) 1_2y E -1 ^ 2 . 


At the lowest energies (ax » v/e~) 


J(E) - BE' 


,-1/2 


So, in the low energy limit,, the new spectrum is a power law with 
spectral index l/2, regardless of the input spectral index. 

Examining the case where all particles are stored for a 
constant time, it is found that in the low energy limit 
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s 


dE 

dt 


dE /— dE „ _ 

•v — a, /e __ <v -C E . 
dx dx 


Using a similar argument to that given above, it is found that for 
an incident power law spectrum Jq(Eq) = A Eq \ at the lowest 
energies the new spectrum is 

J(E) - exp[Ct(l-y) ] E" Y 
i.e., the spectral index is unchanged. 


Model Calculations 

We are now in a position to determine the effects of propaga- 
tion and storage on the spectra of the particles. A number of 
model calculations have been made involving releasing an initial 
spectrum of particles (for simplicity assumed to be a power law 
j(E) = Jq E Y ) at a particular point in the corona, propagating 
them to another location, storing all particles there for a given 
time , and then allowing the particles to propagate outward to about 
10 Rq (the last value was chosen because here the electron density 
is about 5 orders of magnitude lower than the coronal maximum and 
any further propagation of the particles will contribute negligibly 

to the ionization loss). The corona has been assumed to be isothermal, 

6 

with an electron temperature of 1.5 x 10 °K. An empirical fit to 

Q 

the coronal density profile of Allen [1963], n g = 3-6 x 10 J 
exp(l5Rg/R) cm ^ .was- used. Energy loss was determined by numerically 
integrating equation ’ (^-l) along the particle trajectories. 

In figures. U-3 and U-U results of the model calculations are 

—3 

shown for various ions, for. an initial power law spectrum J(E) a. E 
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—2 

and J(E) v E , respectively, propagated outward from an injection 

g 

point at 1.02 R^ to a storage area at 1.4 R^ (where n^ = 1.4 x- 10 
—3 

cm ) , trapped at 1.4 R^ for 1000 seconds, and then propagated 

outward to 10 R . This gives a time-density product of about 
11 -3 

2 x 10 sec cm , which is of the order required to reestablish 

charge states characteristic of 1.5 x 10^ degrees. The particular 

values of storage radius and time do not affect the limiting slopes 

of the curves in figure 4-3, but only the time in which these 

limits are reached; the curves could as easily describe a spectrum 

^ 3 

stored at n = 5 x 10° cm~ (about 2 R_) for 10 hours, 
e ® 

2 

The factor Q /M in equation ( 4—1 ) is apparent in the difference 
in degree of flattening for alphas and heavier nuclei. Clearly a 
"kink" occurs in the spectrum, at the same energy for all species. 

In addition, the calculated energy spectra of both protons and 
helium -show a mild turn-down. 

Comparison with Observations 

We note that the energy at which the bend in the spectra 

in figures 4-3 and 4-4 occurs is about a factor of 10 higher than 

the energy at which the spectra are observed to turn down in the 

events of this study (see figure 3-l4) 

The energy at which the bend occurs, E^ is, to first order, 

a function only of the parameter u in equation (4-l), or of 

1/2 

(E./T )• . In order to make the calculated spectra consistent 

i e 

with observations, E must be lowered by about a factor of 10. 

13 

One way to accomplish this is by lowering the temperature of the 
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storage by a factor of 10. Limits which have been put on the 

temperature of the storage region from the observed charge states 

of the ions, however, are greater than 10^ °K, far too high to 

produce the flattening at the observed energy. Storage at low 

temperatures {with correspondingly low equilibrium charge states) 

followed by stripping may be ruled out by the following argument: 

A turndown in the spectra of CNO which is more pronounced than for 

Q 2 . 

alphas, requires that in the storage region the parameter — m 

equation (4-l) must have been greater for CNO than for alphas, 

assuming all species are stored for the same time. Since for 
D 2 2 

alphas, ; = 1, Q > 12, or .Q > 3.5 for carbon, which implies that 
in order for these high charge states to be maintained, the 
temperature in the storage region must again have been > 8 x 10 °K. 

Therefore, coronal propagation and storage alone could not have 
produced the spectral flattening at the energies where it is observed . 
If the spectral turndown is the result of particle trapping and 
ionization loss, then some additional mechanism must also operate 
before we observe the particles. 

A Model for Particle Propagation 

On the basis of our discussions, it is suggested that the 
observed charge states and spectral bend come about in the 
following manner: Particles (alphas and heavier nuclei) are 

accelerated in a source region to a power law spectrum in kinetic 
energy. The particles are trapped near the source region at 
temperatures T g ^ 1.5 x 10^ °K characteristic of a quiet corona. 

The particles lose energy through ionization loss, their spectrum 
flattens at about 800 keV per nucleon and the ions establish charge 
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state equilibrium with the ambient plasma. From the observations 

of the observed correlation between the degree of flattening and 

the low energy proton flux, it is concluded that either the 

storage time or the storage electron density is proportional to 

the "size" (i.e., total energy) of the flare. The particles are 

then diffusively propagated outward to the earth, losing energy 

by adiabatic deceleration. Here it is assumed that the diffusion 

coefficient is independent of particle parameters [Ma Sung et_ al , 

1975], so that the adiabatic deceleration does not affect the 

shape of the spectrum, but simply moves it to lower energies 

[Venkatarangan and Lanzerotti, 1975]- Under these conditions, the 

required energy losses may be produced by a diffusion coefficient 

20 , 2 — 1 

of the order of 5 x 10 cm sec [Parker, 19&5] which is consistent 
with our observations for at least one event (1974 June 10), and 
with the observations of others [Lin et al, 1974; Lupton and Stone, 
1973] at low energies. Thus, the "bend" in the spectrum of CNO 
is moved down to the energy at which it is observed. Statistical 
errors in the helium observations, unfortunately, prevent an 
unambiguous identification of the predicted bending in the helium 
spectra. 

It should be noted that the value of tc u required above is 
only a average value. A diffusion coefficient which is very small 
close to the sun could have a large effect on energy loss, as most 
of the loss would occur -there [Jokipii and Owens, 1975]* 

Because the theoretical' results derived for alpha particles 
apply equally for protons , ' there should also be a bend in the proton 
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spectrum to match that of the alphas. Since the lowest observed 
energy for protons is even higher than that for alphas, it is not 
expected that this bend would be seen in the data presented here. 

It should also be noted that all species appearing in 
Figure 4-3 were initially equally abundant. After storage, the 
heavy ion abundances are depressed relative to helium. Thus, 
in order to explain events which show enhanced heavy element 
fluxes, it is necessary either to have a "super enriched" source, 
or to invoke some mechanism which preferentially accelerates the 
heavier ions. 

The mechanism of thermal diffusion in the corona, as 
proposed by Nakada [19693, can produce a source region which is 
enriched (to practically any degree desired) in heavier elements, 
depending on thermal gradients in the lower corona. The ionic 
charge states which are observed are consistent with those 
required by his model. For example, the iron charge states 
observed in the May l4 event would produce an enhancement of the 
Fe/O ratio of 'u 10 for a coronal temperature gradient of l°K/cm. 

The Helium/CNO Ratio 

The time profile of the alpha/CNO ratio must now be explained. 
Although the time variations of the ratio during an event are 
consistent with the observations of Van Allen et_ al [1974], our 
interpretation of the cause differs . Van Allen et^ al interpret 
their results as the effect of a rigidity dependent diffusion 
coefficient on partially stripped carbon and oxygen. The results 
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of this study indicate that the rigidities of low energy carbon 
and helium at the same energy per nucleon differ by only 20-30$, 
and that of low energy oxygen is ^ 30$ higher than low energy 
helium. It is unlikely that such a small change in rigidity could 
produce the required large change in the interplanetary diffusion 
coefficient. At the source , however, the composition is a strong 
function of temperature [Nakada, 19&9; see especially his 
figures 5 and 6]. A gradual cooling of the source region 
throughout the event could produce the observed alpha/CNO profile. 
As the source region cools , the source becomes less enriched in 
heavier nuclei (relative to helium), and the alpha/medium ratio 
in the accelerated particle population increases. This source 
region cooling has been observed in soft x-ray emission from 
flares [Craig et al, 1973; Peterson et al, 1973]. Although these 
observations give timescales for cooling of the order of 1 hour, 
they extend only to about 8 x 10^ °K. At lower temperatures the 
cooling time could be longer. We note that this interpretation 
is consistent with the observations and interpretation of Gold 
et^ al [ 1975 ] and Roelof et_ al [1975] of "coronal control" of 
cosmic ray events. 

The same mechanism of thermal diffusion at the source could 
explain the variability of the average alpha/CNO ratio from event 
to event. The heavy ion enrichment depends on thermal gradients 
in the lower corona. As these gradients would be expected to 
vary greatly from flare to flare (due to differences in flare 
geometry), the alpha/medium ratio would also be expected to change 


from event to event. 
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Summary 

The importance of observations of low energy ions from 
solar events as a probe of flare conditions cannot be overemphasized. 
From analysis of ten events, we can conclude: 

1) Low energy carbon and oxygen- appear in ionization states 
which bear the "imprint" of temperatures characteristic of the 
quiet lower corona, and show remarkable consistency from event to 
event . 

2) The very consistency of these charge states implies 

that the low energy ions do not reach us directly, but are trapped 

in regions of electron temperature T g ^ 1. 5 x 10^ °K for some 

. 11 -3 
time - i.e., a time-density product of ^5 x 10 sec cm . 

3) During this trapping the ions spectra are degraded 

by ionization energy loss. This degradation is consistent with 
the observations if strong adiabatic deceleration takes place, 
with the interplanetary diffusion coefficient having the value 
measured in one event. 

b) A time variation of the alpha/medium group ratio has 
been found which is consistent with b.- variable source enrichment 
due to source region cooling. 

5) The ionic charge states of low energy iron have been 
measured for one event , and here too the ionization states are 
characteristic of an electron temperature of 1.5 x 10^ °K. 

Additional event statistics from current experiments, and 
results from even lower energies with greater resolution in future 
experiments, promise a wealth of new data on solar particle 
acceleration and particle propagation in the corona. 
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- APPENDIX A - 

THE UNIVERSITY OF MARYLAND 
ELECTROSTATIC ENERGY-CHARGE ANALYZER EXPERIMENT 

The measurements of the energy and charge spectra of low 
energy interplanetary particles requires two parameter analysis 
from which the particle * s total energy and charge can both he 
determined. 

A technique commonly used consists of measurement of the energy 
loss and total energy for particles using a combination of thin and 
thick detector elements respectively. Telescopes employing the 
so-called "dE/dx vs_ E" method have achieved excellent chemical and 
isotopic resolution by using low noise silicon detectors as the 
sensitive elements of the sensors. In the past the low energy limit 
for two-dimensional dE/dx vs_ E analysis could not be extended 
below several MeV per nucleon because of lack of suitably thin dE/dx 
elements. 

The usefulness of the dE/dx vs_ E analysis technique has been 

extended to energies below 60 keV per nucleon for iron by using a 

thin-window proportional counter as the dE/dx detector element. The 

Ultra Low Energy Telescope (ULET),- included as part of the University 

of Maryland Experiment on the IMP 7 and 8 spacecraft uses a low mass 
2 

(^330 pgm/cm ) total thickness gas proportional counter as the energy 
loss rate, or dE/dx, detector and thus has a low energy limit for two 
dimensional analysis of 300 keV per nucleon for protons and 200 keV 
per nucleon for alpha particles. 
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To extend the range to even lower energies, and to determine the 
particle’s ionic charge requires a different technique. The Electrostatic 
Energy-Charge Analyzer is the main sensor of the University of 
Maryland IMP 8 experiment. This device uses electrostatic deflection 
vs total energy measurements to determine the ionization state and energy 
of particles in the range of 37 to 1250 keV per charge (for charge 
Q _> 2 particles). The sensor 3s similar in principles of operation to 
the instrument flown on IMP 7 [Fan et al, 1971], "but is primarily 
designed to investigate positively charged ions and has a lower energy 
limit — a factor of three "below that of the IMP 7 experiment. In this 
appendix the sensor and the electronics used for signal processing and 
data handling are described. 

The Deflection System 

The deflection d of an ion travelling through an approximately 
transverse electrostatic field generated "by a potential V is related 
to the particle's kinetic energy T and its ionic charge Q through 
[Spiker, 1970] 


d 


= gy S 
(y+1) t 


(A-l) 


where g is a constant dependent on the deflector geometry and y is 
the relativistic factor. Thus, by measuring simultaneously a 
particle's deflection in a known electrostatic field and its total 
kinetic energy one determines its ionic charge Q. 

In the IMP instrument the deflection fields are created by 
applying voltages to two flat non-parallel magnesium plates, mounted 
on stand-off insulators inside a ground cavity (Figure A-l). One 
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plate is maintained at a potential of + it. 5 kV relative to the 
ground plate, the' other at + 3.5 kV. Three field regions are thus 
created in the deflection cavity with potential differences of it. 5 
kV (top), 11 kV (center), and 3.5 kV (bottom), corresponding to 

t 

fields at the front edge of deflection plates of 7*25 kV/cm," 11 kV/cm 
and 1.75 kV/cm. These fields decrease as one moves away from the 
front edge due to the increased separation of the plates. However, 
this does not significantly affect the validity to first order of the 
expression (A-l) above. A passive collimator at the entrance of the 
cavity defines the arrival directions of incoming particles. 

A charged particle entering the sensor is deflected in one of 
these three field regions. Its deflection is determined by the 
location of that one solid state detector in the arrby of eight which 
is triggered as the particle strikes its surface after exiting the 
deflection cavity. Simultaneously, pulse-height analysis of the trig- 
gered detectors output signal determines the particle's total energy, 
differential energy range per charge of each detector is defined by 
its position and width (figure A-2) and does not depend on electronic 
threshold levels. Using eight discrete detectors we obtain eight 
differential energy channel analyzers (seven for positive ions and 
one for electrons) simply counting the number of particles which 
trigger each of the eight detectors. The charge and finer energy 
resolution within each of the 7 positive ion channels is achieved by 
pulse height analysis of the signal's of that detector. 


The 



The Detector Array 


The front of the detector array (Figure A-3) consists of seven 
thin (100 p) silicon surface barrier detectors designated PI to PU 
and P6 to P8 in order of increasing energy of the particles able to 
reach them. In addition, a thick (1000 y) Li drifted Si [Si (Li)] 
detector (called N) is included in the array at a position to inter- 
sect electrons. The signal from the N detector 'is not pulse-height 
analyzed — it serves only as a rate trigger. Detector P5 is covered with 
an aluminum foil just thick enough to stop protons deflected to it. 

This detector is identical to P2 and P6 and its purpose is to measure 
the shape of the background energy spectrum and the background rate. 
Continuous monitoring of the background spectrum and rate are essential 
because the background spectral shape was found to change during solar 
flare events. The detector energy per charge ranges and geometrical 
factors are listed in Table A-I. 

Behind this front layer is an array of six thick Si(Li) detectors, 
positioned so that particles which penetrate one of the front detectors 
must also strike one of these. All back detectors are connected in 
parallel and are referred to as B. The requirement for pulse-height 
analysis of a given P detector P ± is ? i • B, since positive ions 
reaching these P detectors will be stopped in them. High priority for 
analysis is given to heavy ion events in certain detectors. The N 
rate requires N*B coincidence reflecting the high probability that 
electrons reaching N will penetrate and trigger B. The coincidence 
requirement assures identification of electrons unambiguously and 
serves to reduce background. The B anticoincidence for pulse-height 
analysis also helps in reducing background due to penetrating particles. 
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A plastic anticoincidence scintillator surrounds the detectors, 
preventing analysis of high energy penetrating particles and 
further reducing background. In addition, the scintillator 
counting rate is a measure of the integral proton flux above 
about 18 MeV, with a large geometrical factor (about 1100 cm sterad) . 

The Collimator 

The collimator at the entrance of the sensor cavity overcomes 
one fo the major disadvantages of the electrostatic deflection 
method, namely the inherently small geometrical factor. For deflection 
measurements to be meaningful the entrance angles of the entering 
particles must be severely restricted. This entails very strict 
collimation, with the result that such sensors usually have very small 
geometrical factors. We have overcome this problem to some extent 

1 

by using a collimator with a large open area and many slits. 

The 40 mm thick collimator consists of a total of 2J0 thin 

(0.l44 mm) aluminum plates stacked together and rigidly attached to 

either side of a magnesium frame (Figure A-4). Aluminum was chosen 

because of low atomic number and moderate density, yielding a collimation 

2 

thickness of approximately 5-4 gm/cm for particles with skew trajectories. 
The curved frame (inside radius 270 mm) has 150 plates on the convex 
front side and the remaining 120 on the back. Each stack is secured to 
the frame by screws. 

Each aluminum plate has 1521 slotted openings chemically milled 
through its surface (Figure A-5). These slots are 0.254 mm wide and 
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vary in length from 12.37 to 14.25 mm- Proceeding from the outer- 
most to innermost plate the slots become shorter and the row 
spacings narrower. A total of 27 different plate types are used. 
Stacking these plates creates a series of vertically and radially 
converging tunnels through which particles may pass. The upper 
and lower collimator sections have a vertical acceptance angle 
of 6.10° while the center section has a convergence angle of 2.1+7° • 

The horizontal acceptance angle (in the ecliptic plane) is 'v 60°. 
"Focusing" of particles is achieved by the combined vertical and 
horizontal convergence of the tunnels. The. points of convergence 
lie along three separate ar.cs, one for each electric field region, 

188.5 mm behind the collimator. The dispersion of particles on this 
arc represent the geometrical resolution of the focused collimator. 

This resolution depends on slot width, focal length, collimator 
thickness and stacking accuracy. Resolution of this collimator 
ranged from 1.6 to 1.7 ram (FWHM). The geometrical factors, a 
function of both acceptance angle and detector area, were calculated 
for the various detectors using a Monte-Carlo method and are listed 
in Table A-I, above. 

To alleviate small angle scattering of both light and energetic 
particles in the collimator 30$ of the plates were replaced with 
spacers to form "light traps". Three spacers were included with the 
10 plates of each type for .all 27 different types. 

The Electronics 

A simplified block diagram of the electronics is shown in Figure 
A-6. The electronics serves to count the number of particles which 
strike the various solid state detectors and to measure the energy lost 



in the detectors. It also interfaces with the spacecraft, transfers 
the data to spacecraft for telemetry and monitors the performance 
of the instrument. 

Each detector has a set of linear circuits consisting of a 
charge sensitive amplifier, pulse shaping amplifier and at least 
one amplitude discriminator. The output from these circuits is a 
logic pulse from an amplitude discriminator indicating that a pulse 
exceeding threshold has occurred, and a voltage pulse suitable for 
pulse height analysis whose amplitude is a function of the charge 
deposited on a detector. The amplifier noise level is 50 keV 
(FWHM). 

There are two pulse height analyzers (PHA) in the instrument. 

Each PHA contains a delay line, linear gate, pulse height to pulse 
width converter and a gated clock. The maximum PHA channel number, 
determined by the saturation of the amplifiers, is about 1150. The 
spacecraft telemetry limits the analysis rate to one event , con- 
sisting of 2 PHA's, every 1.28 seconds. The two PHA's analyze 
either coincident pulses from ULET detectors or two single pulses from 
EECA detectors. 

Since there are a total of ten detectors whose outputs are 
pulse height analyzed, an increase in the noise level of one 
detector to a point where noise pulses are analyzed will seriously 
reduce the amount of PHA data obtained from other detectors. To 
minimize this possibility the PHA's are multiplexed between ULET 
and EECA so that only ULET pulses are analyzed for a specified 



fixed time period and EECA pulses are analyzed the rest of the time. 
Either sensor can be selected by a ground command to have 75$ (48 
sequences out of 64) of the PHA time. As a further protection 
against individual detector failure, there are threshold rate 
meters monitoring the pulse rates from amplitude discriminators. 
Whenever the pulse rate from a discriminator exceed the threshold- 
rate, pulses from the detector associated with the particular 
discriminator are no longer pulse height analyzed and therefore 
have no adverse affect on the Head times of other, operational 
detectors. All detectors, except the anticoincidence scintillator 
PM tubes, have such rate meters. The rate thresholds are set at 
a. 50 khz for both ULET detectors and the low energy EECA detector. 
All other rate thresholds are at 10 khz. 

The logic circuits generate the event counting rates and 
determine when and which pulses are to be pulse height- analyzed. 

The event counting rates are summarized in Table A-II. The basic 
requirement for positively charged particle counting rates is that 
the front detector produce a pulse and that neither the B detector 
nor the anticoincidence cup (Al) produces a pulse. The basic 
requirement for the electron rate ‘is that both the front detector 
and the back detector have coincident pulses with no pulse from 
the anticoincidence cup. (The nominal requirement to record a 
coincidence or -anticoincidence is that the pulses are within 0.6 
ysec). Most rates are sectored into four quadrants, according to 
aspect angle of the rotating spacecraft with respect to the sun, 
thus providing anisotropy information. The spacecraft spin rate 
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of 22.9 KPM gives a short time scale for these measurements. For 
the University of Maryland experiment 12 rates are available, 8 
sectored into four quadrants each (SR) and 4 unsectored rates (R). 
Since a total of 18 rates had to he telemetered some rates are 
multiplexed as shown in Table A-II. 

The PHA requirements for the EECA detector analysis are the 
same as for positive ion rates with the additional requirements 
that the analyzer is not busy, the data from registers are not 
being telemetered, no other detector using the same PHA has a 
pulse, the PHA is in the EECA instead of ULET mode, the threshold 
rate from the detector has not been exceeded, and that no high 
priority event is stored in the register. High priority is assigned 
to analyses of charge Q >_ 4 from P2 detector in one PHA and to 
analyses of charge Q > 2 from P7 in the other PHA. When these 
analyses occur a flag is set and the PHA is inhibited from analyzing 
other events until the data is telementered. If the priority flag 
is not set, pulses are analyzed until a. 2-5 msec before telemetry 
and the result of last analysis is transmitted. The priority scheme 
was implemented to obtain PHA data for heavier particles during times 
of high fluxes. With each PHA event a code identifying which 
detector produced the pulse and the angular position of spacecraft, 
to the nearest 45°, at the time the event occurred is also recorded. 

To verify the proper operation and to detect any significant 
drifts in gain or thresholds an in-flight calibrate system is used. 

For a period of a. 82 seconds every a, 3 hours the instrument is 
commanded by the spacecraft into a calibrate mode. During this 
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time pulses are applied to the test inputs of each of the charge 
sensitive amplifiers in turn. The amplitude of these pulses is 
stepped through a predetermined sequence. For each linear channel 
pulses of appropriate amplitude are generated so that the overall 
gain from charge input to channel number output is checked at at 
least two points. To avoid particle data from interfering with 
calibration the PHA input is enabled only for 'u 1+0 usee and the 
calibrate input pulse is generated during this time. 

Analog voltages which monitor performance of some subsystems 
are also periodically read out. The temperature at the EECA detector 
■cup and the detector bias voltage are measured once every 'v 82 seconds. 
The + 3.5 kV supply, the lb. 5 kV power supply, ULET proportional 
counter pressure and the temperature of the thermal valve controlling 
ULET pressure are monitored every 161+ seconds. 

Extensive calibrations were performed on the EECA, using both 
sources and accelerator beams. The first calibrations and design veri- 
fication have been reported by Spiker [ 1970 ]. Additional source cali- 
brations and resolution measurements were reported by Fan et al [l97l] - 
Electron, proton, helium, and nitrogen beams were used to calibrate the 
IMP -8 experiment, both at the University of Maryland and G-oddard Space 
Flight Center accelerators. 
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INSTRUMENT RESPONSE 


In order to make full use of the PHA histograms, it is 
necessary to understand how the instrument responds to a incident 
spectrum of particles, and to deconvolve the effects of dete’ctor 
and electronic noise, deflection dispersion, electronic thresholds, 
energy defect, etc. 

We first define: 

a) e(E) = detection efficiency: the probability that 

a particle incident with energy E strikes 
the detector active area. 

We note that a particle of incident energy 


E does not, in general, deposit its full energy 
in the detector, because of non-ionizing nuclear 
collisions [Campbell and Lin, 1973; Linhard et al, 1963]. 
The deposited energy E’ is related to E by 
E» = fi(Z,A,E) • E 


where 6 is the nuclear defect and 5 _< 1. An 
empirical formula for 5 has been derived and used 
in this analysis : 

<5 (Z,A,E) = 1 - {[l-5.12xlO~ 3 (Z+lO)sin(|^- tt)] x 


.03737 


/ 


Z&nZ 

E 


with E measured in MeV. 
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b) D(E' ,h) dh = probability that deposition of an energy 
E 1 in the detector gives a pulse height 
between h and h+dh. This is ideally a 
Gaussian function of the electronics noise with 


V 


•D(E*,h) = — 


-(E*-h) 2 /2a| 


S2ts a. 


E 


c) P (h) = probability that the detector signal h 

Xt 

triggers the threshold discriminator, allowing 

pulse height analysis to take place. 

= , where t is the threshold and P(x) is 

a E 

the normal probability integral 

p(x) - -i- r .-*!* 




L 


du 


We can now define a detector response function [Berger et al , 1 969 ] 


R(E,h)dh = e(E) D(6*E,h) P t (h) dh 

= probability that an incident particle with 
energy E is detected and gives rise to an 
analyzed pulse height between ’h and h+dh. 

If the incident flux of particles is j(E) (particles/cm -sr- 
sec-MeV), we then obtain a pulse height distribution H(h) 
(particles/sec-MeV) , * • 


H(h) dh = G 


J .5(E) R(E,h) dE • dh 
0 


where G = geometrical factor (cm sterad) 

H(h) dh = # of particles detected/second with pulse 
height between h and h+dh 
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= G f j(E) e(E) D(6-E,h) P (h) dE • dh 

Jo 


Going to the discrete case, for pulse height channel n we define 


H = 
n h 


i A +1 

n+l \i «' h 
n 


H(h) dh 


= particles/sec-MeV in channel n with and h n+1 the lower 
and upper energy limits of channel n, respectively. 


H = - 
n h 


i r n+1 r 

““r — / dh 0 I j(E) e(E) D(E6 ,h) P (h) dE 


n+l n ^ h 


n 


h 


" /-n+i 

-r— / dE j(E) e(E) J — 

~K Jr, v, To 


n+l n ^ o 


E 


— , h 
m=l n+l 


-(ES-h) 2 /20 2 

E P(^) dh 

/ m+1 rn+I -(E6-h) 2 /2<? 2 

dE j ( e)e(E) / — e E P(^) dh 

J h v^2tT o 

£i n E 

m 


E 


h fa o_ 
n E 


E 


N 

or H * V' - 
n h 


m=l 


■ “ h 
n+l n 


E 

/ m+l E6-h 

dE j(E) e(E) [P(— ~) - P( 

E E 

m 


E6-h 


n+l 


•)] 


E 


h -t 

•p(-g — : i 

E 


Recall that e(E) = detection efficiency for incident energy E, 

or the probability that a particle with incident energy E will 

* 

strike the detector. For a perfect electrostatic deflection sensor, 
e(E) should be a rectangular window: 




- 50 ' 


where E . and E are determined hy detector size and position, 
man max * 


and particle deflection y: 



However, in a real sensor, particles of energy E q are 


deflected with some finite dispersion about the mean position 



and thus the detection efficiency looks like this: 


Numerically 



where y Q = and ct d = (.^328l + .00081 y 2 ) 1/2 mm [ Fan et aL, 
o 

197l] is the deflection dispersion. 

Thus ' 

s(e ) = . p(VijaS) 

0 °D °D 
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rE y min 


kg _ 

) r ( E ~ ^max ^ 
°D 


So the pulse height histogram becomes 


_-!X; P r m+l ES-h E6-h 

‘n • £ r^r i iWiV - 1 


m=l n+1 n E 


kg. kQ 

[p(£- Iika. ) . P ( S . P( _sL_) ^ 

°D °D °E 


particles/(sec MeV) in channel n 


If we assume that over the energy interval (E^, j(E) is a 

constant (a good approximation for small channel widths E m+1 - E ffi ) 


_ -i r m+x 

= E ”_-E J 
m+l m E 
m 


(B-l) 


then. 


.m+l E6-h 


E6-h 


H = r — ) 
n h , , -h £ 
n+1 n m 


C i m f w-sr^ - 

1=1 m 


— kg 

Cf( e " ” ln ) - p(S-lla5£)] p(-al) ajj 

D E E 


or H = > R j 
n £—i nm m 
m=l 


/ m+l ES-h j 

[ P (_n) _ P (. 

0 V 


nm h ,,-h 
n+1 n 


S-h ^ - y 

S±l)] [p( E mis.) 

°E °D 


M_ h 

- P<^ SS)] P(-S^) dE 

0 D a E 


(B— 2 ) 
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H 


Thus 



(B-3) 


where H and j are N-dimensional column vectors and E is an NxN 
matrix. 
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- APPENDIX C - 


DETERMINATION OF CHARGE STATE ABUNDANCES 


Suppose that we have a pulse height histogram C which we 

"believe may he represented as a linear superposition -of a set of 

histograms P., where P. is the result of operating on an incident 
0 J 

spectrum of particles of type j with the instrument response function 
B(B-3). As an example, assume that for a particular event a given 
detector's histogram is a composite of the charge states 5,6,7, and 


8. Thus, for a given element (channel) of the histogram C^, 
C i “ a 5 P 5i + P 6i + a 7 P 7i + a 8 P 8i * 


(C-l) 


Our problem is now to determine the .most likely . 

Since there will he statistical errors in the C^, the equality 
in equation (C-l) will not, in general, hold. Define a difference 
measure for channel i 


A. 2 = 

l 


(C. - V a. P ..)* 
1 hi o a 1 


and a difference sum 


v 2 ■ 


By the principle of maximum likelihood, hy minimizing this sum, 
we determine the most likely values of the a, . 

We now perform a standard least squares analysis. At the 
minimum of S, its variation <5S is stationary with respect to changes 
in the a„ 


SS = 0 
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{ -X>. p --) 2 > * o 

i _ 2 a. *r / ,1 ,n 


K i ct. 

l 


thus. 


IV -L 

Et ,! «V-E « 1 V 5 - 


P TT- C - P~ 

i a/ i J J u. 


interchanging the order of summation, 

o 1 »i 


defining 


we have 


K- c - 

1 a. 
i 


P P 
J i a / 

i 


^ \ d a j 


6 = A a 


[which is' a polynomial equation in the a which may he solved by matrix 

* ” . J 

techniques] 
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a=e| 


with e = A ; hence 


a. = E (e^ e K ) = E U iV - E 


p v . c. 

Kx 1 ■ 


jKr 2 

i a. 

i 


(C-2) 


The errors in the a. are given by [Bevington, 1969 ] c = e . 

J B - J J 

tJ 

Unfortunately, the histograms P are functions of the spectral 

J 

index of the incident spectrum (see Appendix B), so the solutions 

of (C-2) are not well defined. However, from rate data it is 

usually possible to estimate the spectral index and use this estimate to fit 

histograms. In addition, it has been determined that in most cases 

there is only a weak dependence of the calculated abundances on spectral 

index used, and it is these ratios that we are usually interested in. 

A computer program, named M3PCT, has been developed which 
performs all necessary operations on the pulse height data. This 
program uses as its input data pulse height histograms punched on 
80 column computer cards. (These cards are punched by another 
program, which is responsible for accumulating the PHA data for 
specified time periods from the experiment encyclopedia tapes . ) 

MSPCT calls subroutine BSUB to normalize and subtract a quiet time 
"background" histogram from the histogram of the time period of 
interest. The background histogram may be normalized in one of 
several ways. Usually, the histogram is normalized to absolute flux 


levels. During time periods when the flux of high-energy protons is 


high the background counting rate changes, and it has been found 
that this "absolute flux" method does not give very reliable results. 
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In these eases, it is possible to specify a certain energy range 

over which no "real" events are expected. BSUB then calculates 

the normalization constant which will minimize excess counts over 

background in this range. The third method is used when it is 

believed that the background spectrum in a detector has some simple 

form — a power law in energy, for example. The measured histogram 

then consists of this power law background plus, over some energy 

range, a number of "real" events. By fitting a polynomial to the 

histogram exterior to those areas where genuine events are expected, 

the resultant fitted background curve may be used to correct the 

observed histogram. The program then reads in the appropriate 

instrument response matrices from magnetic tape (these matrices are 

generated by a subroutine named, logically enough, MATRIX). The 

routine calls a subroutine which generates an "incident" histogram 

of the form (B-l) for the input spectral index, and performs the 

matrix multiplication of this histogram and the response matrix to 

generate the predicted instrument response histogram P . The matrix 

J 

A is generated and inverted, and the a. are computed. Both the data 

J 

histogram (before and after background correction) and the computed 

’ * 

superposition of P. may be plotted. Finally, charge state abundance 
3 

ratios and a "goodness of fit" measure are computed. 

In addition, the routine has the ability to generate and plot 
pulse height histograms which represent "predicted" instrument 
response to an input mix of charge state abundances and spectral indices. 
The dashed curves of Figure 3-2 were generated in this manner. 
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• Table 3-1 

Characteristics of 10 Solar Particle Events 


Event 

Date 

PHA Data 

Time 

Flare 

Importance 

Approx . 
Location 

73/11/4 

11/4 :0000 to 11/5 :1159 

11/3 :0012 

2N 

S18W85 

TV 1/25 

1/25:1200 to 

1/27:0300 

1/22:1915 

IN 

N09W60 

74/5/14 

. 5/14:0700 to 

5/15:0500 • 

. 5/13:2118 

2N 

S13W66 

74/5/31 

5/31:0000 to 

5/31:1200 

5/29:1200 

IB 

S15W63 

74/6/10 

6/9 :0000 to 

6/11:1200 

6/8 :l450 

3N 

S13W64 

74/6/22 

6/22:0000 to 

6/24:2400 

6/21:1822 

IF 

S20W27 

74/6/26 

6/26:0000 to 

6/ 29:1200 

6/24:1515 

IN 

S13W69 

74/7/4 

7/5 :1200 to 

! 

7/8 :2400 

7/4 :1338 

2B 

S16W03 

74/10/14 

: 10/12:1200 to 10/15:2400 

10/11:1010 

IB 

N13W03 

74/11/8 

11/8 :0000 to 11/10:2400 

11/6 :0220 

IB 

_ S12W75 





Table 3 -I I 


Ionization State Abundances 


Date 

C +5 /C 

c +< 7c 

o +6 /o 

0 

0 

0 +8 /0 

Fe +10 /Fe 

Fe +11 /Fe 

Fe +12 /Fe 

Fe +13 /Fe 

73/11/4 

0.25 i.ou 

0.74 -.Ok 

0.72 i.05 

0.20 i.05 

0.07 i.03 





74/1/25 

0.31 i.06 

0.68 i.07 

0.69 i.o4 

0.16 i.06 

0.05 i.06 





74/5/14 

0.23 -.03 

0.76 -.05 

0.66 i.05 

0.33 i.04 

0.00 i.03 

o.i4 i.o4 

0.30 i.06 

0.34 i.07 

0.20 i.05 

74/5/31 

0.20 -.05 

0.79 -.06 

0.70 i.07 

0.23 i-05 

0.06 i.04 





74/6/10 

0.17 -.02 

0.82 i.05 

0.81 i.07 

0.16 i.06 

0.02 i.01 





74/6/22 

0.34 i.07 

0.65 i.02 

0.73 i.06 

0.20 i.03 

0.06 i.05 




/ 

74/6/26 

0.18 -.02 

0.81 i.07 . 

0.67 i.o4 

0.21 i.06 

0.08 i.04 





74/7/4 

0.28 i.06 

0.71 i.03 

0.68 i.05 

0.27 i.o4 

0.04 i.06 





74/10/14 

0.26 i.05 

0.73 i.o4 

0.77 i.03 

0.21 i.05 

0.01 i.02 





74/11/8 

0.22 i.04 

0.77 i.06 

0.71 i.o4 

0.15 i.06 

0.13 i.08 





All Events 

0 

« 

+ 1 

CM 

. 

O 

0.75 i.05' 

0.71 i.05 

0.22 i.05 

0.05 i.o4 

o.i4 i.04 

0.30 i.06 

0.34 i.07 

0.20 i.05 

Solar Wind* 

3.) 


0.77 

£0.15 


0.24 

0.21 

0.14 


Solar . v 
Corona ' 

0.186 

0.813 

0.759 

0.234 

0.008 

o.i4i 

0.257 

0.316 

0.155 


a) Same et_ al [1968] 

b) Jordan [1969] for T e = 1*58 x 10 6 K 



Table 3-III 


Composition and Spectrum Data for Ten Solar Particle Events 

Event Averaged Values 



Relative Abundance 

Alpha Spectral Index 

CEO Spectral Index 

proton flux 

ULET 

Abundances 

Date 

a/CNO 

Y cc 

y cno 

(cm^-sr-sec-MeV) b 

(600 keV/nuc) 


(40 keV/nuc) 

40-?80 keV/nuc 

20r40 keV/nuc 

130-220 keV 

C/0 

Fe/0 

73/11/4 

14 

1.9 - 0.2 

1.1 - 0.4 

340 

0.43 

<0.01 

74/1/25 

17 

2.4 t 0.3 

0.9 - 0.5 

1000 

1.3 

0.14 

74/5/14 

2 

O 

+ 1 

CO 

CVJ 

2.4 1 0.3 

65 

0.39 

2.5 a) b) 

74/5/31 

17 ' . 

2.6 i 0.3 

1.5 ^0.4 

100 

0.27 

0.09 

74/6/10 

3 

0.4 t 0.5 b ^ 

0.15^0.5 

1200 

0.46 

0.20 

74/6/22 

7 

2.7 i 0.4 

2.0 i 0.5 

200 

0.64 

< 0.01 

74/6/26 

l4 

3.1 - 0.5 

1.1 i 0.5 

' 250 

0.79 

0.16 

74/7/4 

16 

1.8 - 0.3 

-0.2 ~ 0.5 

12000 

0 . 4 o 

0.36 

74/10/14 

20 

3.3 - 0.4 

1.2 t 0.5 

230 

0.31 

0.17 

74/11/8 

8 

1.7 - 0.2 

- 0.5 - 0.75 

2500 

0.46 

0.05 


a) "iron rich" [Gloeckler et al, 1975b] 

b) alpha turndown observed in this event 
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TABLE A-l 

Detector Energy-Charge Ranges 


Detector 3- Charge 

Threshold 


Energy 

Range 


Geometrical 
Factor (cm^ sr) 


PI 

Q 

> 6 

37-55 keV/Q 

.02 . 

P2 

Q 

>_ 2 

65-100 keV/Q 

.03 

P2' b 

Q 

>_ 4 

65-100 keV/Q 

.03 

P3 

Q 

> 1 

130-210 keV/Q 

.03 

p4 

Q 

_> 1 

160-230 keV/Q 

.02 

P5 

Q' 

> 1 

background > l60 keY 

— 

P6 

Q 

>_ 1 

390-600 keV/Q 

.018 

p6 ,c 

Q 

4 

390-600 keV/Q 

.018 

P7 

Q 

> 1 

550-900 keV/Q 

.03 

p j r ^ 

Q 

1 2 

550-900 keV/Q 

.03 

P8 

Q 

> 1 

740-1200 keV/Q 

.018 

N-B 

electrons 

600-860 keV 

.018 

A1 

, protons 

> 18 MeV 

1100. 

A2 

protons 

> 6 MeV 

675. 


discriminator rates for all detectors are generated and read out 
through telemetry. In addition, all detector signals except N*B 
and A1 and A2 are pulse height analyzed. 


Denotes pulse height analysis priority event . 


c 

For detectors listed twice the prime indicates a higher threshold 
for rate or pulse height analysis. 


TABLE A-II 
Event Counting Rates 


Rate 

Multiplex 

Step 

Time in 
Mux Step (sec) 

Logical Equation 

Accumulation time/ 
readout ( sec ) . 

Period between 
readouts (sec) 

R1 



p4 .( b+brd )-ai 

5.12 

5.12 

R2 



P7* • (B+BRD) *A1 

5.12 

5.12 

R3 



P2' •{ B+BRD )-Al 

5.12 

5-12 

R4 

0 

61.44 

A1 . 

10.24 

10.24 


1 

20.48 

,A2 

10.24 

10.24 


2 

61.44 

B 

10.24 

10.24 


3 

20.48 

A2 

10.24 

10.24 

SRI 

0 

61.44 

No rate ( sync ) 

— 

20.48 


1 

20.48 

P5* ( B+BRD )-Al 

4.59/sector 

20.48 


2 

81.92 

P7* (B+BRD) -A1 

4.59/sector 

20.48 


3 

81.92 

P6 • ( B+BRD ) • A3. 

4.59/sector 

20.48 


4 

81.92 

P3- (B+BRD) -A1 

4.59/sector 

20.48 

SR2 



P6'- (B+BRD) -A1 

4.59/sector 

20.48 

SR3 


. 

N- (B+BRD) -A1 

4.59/sector 

20.48 

SR4- 

0 

61.44 

D2-A2 

4.59/sector 

20.48 


1 

20.48 

D1-A2 

4.59/sector 

20.48' 

SR5 



P8(B+BRD) -A1 

4.59/sector 

20.48 

SR6 



D1.D2-A2 

4.59/sector 

20.48 

SR7 



P2( B+BRD )A1 

4 . 59/sector 

20.48 

SR8 



PI ( B+BRD )A1 

4 . 59/sector 

20.48 


NOTE: BED true when B threshold rate is exceeded. 
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figure CAPTIONS 


Figure 2-1 


Fi glare 2-2 


Figure 2-3 


Figure 3-1 


Figure 3-2 


Ionic charge state resolution of IMP-8 detector 
P3. This detector nominally responds to ions 
with energies in the range 130 to 220 keV per 
charge . Note that , while the energy hands, for 
ions of charge 5, 6, 7> and 8 overlap, there is 
a range where only ions of charge 5 appear. By 
means of a simple iteration procedure, it is 
possible to calculate abundances for each ionic 
charge state. 

Cross section of the Max-Plan ck-Insti tut* s Ultra 
Low Energy Telescope (ULET) as flown on IMP-8. 

The AE element is a thin-vindow flow through 
proportional counter. 

The IMP-8 spacecraft orbit in the solar ecliptic 
(SE) coordinate system, plotted from ephemeris 
tapes for early May, 197 1 *- The upper p anel shows 
the distance from the earth-sun line 
vs_ distance from the earth (in earth radii) SE 
along the earth— sun line (X _) . The nominal position 
of the earth's bowshock and magnetopause are shown. 
The lover left-hand panel shows the projection 
of the orbit into the ecliptic (X -Y ) plane. 

The sun-earth-vehicle (SEV) angle is-indicated. 

The lower right-hand panel shows the projection 
in the Y gE ~Z plane, normal to the ecliptic. 

Courtesy of Dr. J.D. Sullivan, MIT, 

a) Uncorrected pulse-height histogram from detector P3. 
The smoothed, measured background distribution is 
shown as a dashed line. 

b) Pulse-height histogram corrected for background. 

Note -the characteristic triple-peaked form of the 
distribution. The rightmost peak corresponds to 
carbon and oxygen in the charge group 5-8, or iron 
with 8 <_ Q <_ 12 . 

a) Corrected pulse-height histogram from detector P3 

for the- "iron rich" event of Gloeckler et al [l975b] . 
The dashed curve represents a computed fit of J 

ionic charges 8 to 12 of iron and 5, 6, 7, and 8 of 
carbon and oxygen to the observed distribution. 

b) The same histogram + ag a). Here it has been assumed 
that all iron is Fe . The required Fe peak 
(dashed curve) is clearly not present. 
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Figure 3-3 


Figure 3-4 


Figure 3-5 


Figure 3-6 


Measured iron and helium fluxes from ULET and the 
EECA during the 1974 May l4 event . Iron has been 
assumed to have ionic charges 8-12 in the EECA. 

The agreement between the two experiments is very good. 

Fluxes of 40, 80, and 150 keV per charge heavy ions vs 
time in the 1974 July 4 event. Note that the fluxes 
' from the two lowest energy channels overlap , 
indicating a strong bend-down in the energy spectrum. 

The he lium/ medium-group' ratio at 80 keV per charge 
plotted as a function of time, for the July 4 event. 
During both large particle increases in Figure 3-4, 
the ratio starts low and increases throughout the event. 

- The event-averaged energy spectrum of. helium during the 
1974 July 4 event. 


Figure 3-7 


Figure 3-8 


Figure 3-9 


Figure 3-10 


Figure 3-11 


Figure 3-12 


Figure 3-13 


Figure 3-14 


The event-averaged energy spectrum for the charge 5-8 
group during the 1974 July 4 event. Note the spectral 
flattening below 100 keV per nucleon. 

The mean ionie charge states of carbon and oxygen as a 
function of energy for the 1974 July 4 event. There 
is no statistically significant variation with energy. 

Proton, helium, and medium-group differential energy 
spectra for the 1974 June 10 event. In this event, 
both the helium and the medium-group spectra show a 
flattening at lower energies . 

The mean ionic charge states of carbon and oxygen as a 
function of energy for the 1974 June 10 event. Again, 
there is no significant energy dependence. 

Time variation of the spectral index y (j ^ E~ Y ) for 
heavy ions between 40 and 80 keV per charge during the 
1974 June 10 event., The spectra^ index, after velocity 
dispersion effects no longer are important, remains 
at a low, more-or-less constant value. 

The helium/medium-group ratio at 80 keV per charge, for 

‘ the June 10 event. '-The trend is for the ratio to 
increase with time during the event . 

Proton, helium, and medium-group differential energy 
spectra for- the 1974 June 22 event. This event shows 
no spectral flattening (compare to Figures 3-9 and 3-14). 

Proton, helium, and medium-group differential energy 
spectra for the 1974 November 8 event. Here a strong 
bend-down appears in the heavier nuclei, but there is’ 
none apparent in the helium spectrum. 
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Figure 3-15 


Figure 4-1 


Figure -4-2 


Figure 4-3 


Figure 4-4 


Figure A-l 


Fi glare A-2 


Figure A-3 


Figure A- 4 


Medium-group "event-averaged’' spectral index (40-80 
keV per charge) vs_ "event-averaged" proton flux 
(130—220 keV). Note the strong correlation of 
(-spectral index) with, event size. 


'The expected charge states of carton, oxygen, and iron 
as a function of temperature [Jordan, 1969] , with 
the experimental results shown as a dashed line . For 
the theoretical points , a tar represents<the tempera- 
ture range over which more than half of the ions 
have that ionization state; a point represents the 
temperature at which that species is most abundant. 
After Uakada [1969] . 


The relative ionization loss rate term 
plotted vs u = T9- 5 


I E. 

1 

T 


erf ( 


U 1 in ^exp(~u 2 
- ±.ld u 


The function has the expected 


form for ion velocities near the plasma electron velo- 
city. 


Results of model ionization loss calculations for 
various ions . All species were initially equally 
abundant , with a power law spectrum J E ^ . The 
particles were released at 1.02 R , stored at 1.4 R 

_ o ® 0 

( n e ~ 10 8 cm ) for 1000 seconds, and then propagated 
to 10 Rq. The corona .has teen assumed to be isothermal, 
with an electron temperature of 1.5 x 10 6 °K. 

The same as^Figure 4-3, except that the initial spectrum 
was J a, E 2 . The bending occurs at the same energy in 
both cases . 


Cross section of the University of Maiyland Electrostatic 
Energy-Charge Analyzer, as flown on the IMP-8 spacecraft. 

Energy ranges of the IMP-8 experiment detectors for 
various’ incident ionic charges. 

(Upper panel) Front view of the IMP-8 detector array. 

This array consists of eight 100 y totally depleted 
surface barrier detectors and one 1000 y Si (Li) 
detector (N). 

(Lower panel) Back view of the IMP-8 detector array. 

Here groups B1 and B2 are shown. Electrically, al 1 
of these detectors are wired in parallel and collectively 
referred to as B. 

The collimator assembly. The stack consists of 270 thin 
aluminum plates, mounted to a curved magnesium frame. 
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Figure A-5 


Figure A -6 


The innermost and outermost collimator plates . 

Proceding from the outermost to the innermost 
plate the slots become shorter and the row spacing 
smaller. This creates a series of vertically converging 
radial tunnels when the plates are assembled to the 
frame . 

A simplified block diagram of the EECA instrument as 
flown on IMP-8. All linear and logic circuits are 
shown, which are required for pulse -height analysis, 
rate generation, and internal calibration. 
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FIGURE 2-1 
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FIGURE 2-2 
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figure 3-3 
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FIGURE 3-7 
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TABLE 3 


Effect of Antifoams on Foam Height^ 
With. Ivory and Palmeto Based Wash Waters 



Foam Height at Various Dosages (cm) 


200 ppm 

300 

ppm 

400 

ppm 

500 

PS 

kJ 

3 


15 sec 

60 sec 

15 sec 

60 sec 

15 sec 

60 sec 

15 sec 

60 sec 

Palmeto Soap Based Wash Water 






AF-60 

0.1 

0.1 

m 

0.1 

- 

- 


- 


0.2 

0.1 

m 

0.0 

“ 

— 

“ * 

— 

AF-72 ( 1 2 3 ) 

0.1 

0.1 

0.0 

0.0 

- 

- 

- 

- 


0.1 

0.1 

0.0 

0.0 





AF -93 (2) 

0.1 

0.1 

K« 

mm 

- 

- 

- 

- 


0.1 

0.1 

mm 

mm 

— 

*" - 



Ivory Soap Based Wash Water 




■I 

m 

(3) 


— 

0.5 

0.2 

— 

— 



SWS-211 

0.6 

m 

0.5 

0.2 

- 

- 

SB 

E9 




0.4 

0.2 

0.4 

0.2 ‘ 

0.5 

0.2 

AF-72 

- 

- 

0.4 

0.2 

0.5 

0.2 

0.3 

- 




0.5 

0.2 




_ 

AF-93 










- 

- 

0.5 

0.2 


- 






0.5 

0.3 


_ 


_ 

AF-60 










” 

** 

0.4 

0.3 


~ * 




(1) By NASA "shake test" 

(2) General Electric 

(3) SWS Silicones 
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TABLE 4 


Effect of Antifoams on Dynamic Foam Height 
Using Miranol JEM Based Wash Water 


Antifoam 

Foam Height (cm) Vs. Time Resulting 
From Various Dosages of Antifoam 

50 ppm 

200 ppm 

1 min. 

2 min. 

5 min. 

1 min. 

2' min. 

5 min. 

AF-75 (1) 

pfiiB 

1/2 

1/2 

1/2 

1 



- 

- 

- 

DB- 110 

1/2 

1/2 

1/2 

1/2 

DB 

- 

- 

- 

(3) 

Colloid 680 

- 

- 

- 

0 

1/2 

1/2 

1/2 

1/2 

1 

Control, no 
antifoam 

4 

7 

13 

- 

- 

- 


(!) General Electric Company 

(2) Dow Corning Corporation 

(3) Colloids, Inc. 
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TABLE 5 


Effect of Antifoams on Dynamic Foam Height 
Using Palmeto Based Wash Water 


Antifoam 

Foam Height (cm) Vs. Time Resulting 
From Various Dosages of Antifoam 

200 ppm 

300 ppm 



1 min. 

2 min. 

5 min. 

1 min. 

2 min. 

5 min. 

AF-75 

TOfU 

2 

2 

•3 

4 

1/2 

0 

1 

1/2 

EB 

(1) 

1-1/2 

2 

2 

1/2 

1/2 

1/2 

SWS-211 

1 

1 

1-1/2 

1/2 

1 

1 

DB-110 

2 

2 

4 

1 

1-1/2 

mm 

1 

1-1/2 

4 

1/2 

1 

Wm 

Control, no 
antifoam 

48 cm (5 minutes) 





(!) SWS Silicones, Division of Stauffer Chemical 


































TABLE 6 


Effect of Antifoams on Dynamic Foam Height 
Using ML-11 Based Wash Water 



Foam Height (cm) Vs. Time Resulting From Various Dosages of Antifoam 

Anti foam 

100 ppm 

150 ppm 

200 ppm 

300 ppm 

1 min. 

2 min. 

5 mim 

1 min. 

2 min. ' 


1 min. 

2 min. 

5 min. 

1 min. 

2 min. 

5 min. 

(1) 

AF-93 

MM 


4 

3 

- 

- 

- 

0 

0 


1 

1-1/2 

0 

0 

1 

SWS-211 

- 

- 

- 

- 

- 

- 

1 

Ilj 

2 

0 

0 

1 


- 

- 

- 


- 

- 

1-1/2 

1 

2-1/2 

- 


- 


l 

2 

5-1/2 




0 

1 

1-1/2 

0 

0 

1 

AF-75 





- 

- 

0 

1 

2 

- 

- 

- 

Ferric 



■* 

1 

1-1/2 

1-1/2 







Chloride 













Control, no 




15 

25 

43 







antifoam 














(1) General Electric Company 































TABLE 7 


Effect of Antifoams on Dynamic Foam Height 
Using Ivory Soap Based Wash Water 


Antifoam 


Foam Height (cm) Vs. Time Resulting 
From Various Dosages of Antifoam 


100 ppm 

200 ppm 


300 ppm 



1 min. 

2 min. 

5 min. 

1 min. 

2 min. 

5 min. 

1 min. 

2 min. 

5 min. 

AF -93 

0 


2 

0 

1/2 

1/2 

0 

1/2 

; . - ' 

0 

H9 

2 

1 

1-1/2 

3 

- 

- 

| 

AF-75 

2 

2-1/2 

3 

0 

1/2 

1 

Ht 

- 

- 

2 

3 

5 

1 

1 

1-1/2 

- 

- 

- 

AF-72 (1) 

1/2 

1 

2 

0 

0 

1/2 

0 

0 

0 

1-1/2 

2 

3 

1 

1 

1-1/2 

1 

1 

1-1/2 

DB-110 

1 

2 

5-1/2 

0 

1/2 

1-1/2 

1/2 ‘ 

1/2 

1-1/2 

1 

2 

5-1/2 

1 

1 

2-1/2 



- 

Ferric 

Chloride 

- . 

- 

- 

0 

- 

0 

- 

- 

- 

Control, no 
antifoam 

27 

37 

48 

- 

- 

- 

- 

- 

- 


(1) General Electric Company 






























III. EVALUATION OF REMOVAL AGENTS AND PROCESSES 


Under Contract NAS 9-13536 - Wash Water Solids Removal System Study - 

we demonstrated the feasibility of chemical precipitation using ferric chlor- 
ide. During the follow-up program under Contract NAS 9-14518, we continued 
our evaluation of removal agents using synthetic wash waters based on several 
soaps including Ivory and Palme to . Under the present program we continued to 
screen removal agents for various wash water/soap combinations. 

We have continued our precipitation work using the traditional "jar ex- 
periments". General procedures have been described in "Wash Water Solids 
Removal System Study", Final Report, under Contract NAS 9-13536, July 1974. 
Coagulation work was done with a Hach Floe Tester, essentially a gang stirrer, 
which consists of six magnetic stirrers connected in tandem and set to oper- 
ate at 60 rpm. Coagulants were added with a dasher mixer; mixing time was 
30 seconds. The slow mix time was 15 minutes at 60 rpm. 

PALMETO BASED WASH WATER 

A synthetic wash water was prepared as described on page 4 using- 1000 
ppm (soap solids) of Palme to. Preliminary precipitation experiments using 
aluminum sulfate were not promising (Table 8) . There was very little coagu- 
lation of the soap in the range 200-370 ppm aluminum sulfate. We evaluated 
lower dosages of aluminum sulfate (160-190 ppm - Table 9) , but again the 
results were not promising. No further work was done with this combination 
of soap and coagulant. 

NEUTROGENA BASED WASH WATER 

During earlier work under Contract NAS 9-13536, we found that 160-190 
ppm of ferric chloride would generally remove 60-70 percent of the soap. The 
removal estimate was based, however, on total dissolved solids, which later 
in the program was suspected of being somewhat inaccurate. To better define 
what is happening during the Neutrogena coagulation, we refined our analysis. 
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According to U. S. Patent 2,820,768 and the McDonnell Douglas Report 
covering OSW Contract 14-30-3062, Neutrogena Bar Soap is composed of the 
following : 

1. Sodium soap (saponified tallow, coconut oil, and 

castor oil) — 35—40 percent 

2. Triethanol ammonium salt of stearic acid - 35-40 percent 

3 . Triethanolamine ) 

. ) Relative amounts not 

) specified by manufacturer 

5 . Perfume ) 

Based on chemical composition, this represents approximately 55-60 per- 
cent minimum quantity of fatty acid (normally mixtures of C-14 to C-18 fatty 
acids such as stearic and palmitic) . 

Our chloroform/IR analysis indicated approximately 50 percent fatty 
acid in Meutrogena. This figure is perhaps low, but not unreasonable, as- 
suming that the Neutrogena Corporation is operating within the claims of 
their patent. 

We prepared fresh wash water with 1000 ppm soap. Since there had been 
no evidence that the electrolytes in the synthetic wash water were playing 
a part in the coagulation , we left them out in this case. This was done to 
facilitate the analysis. 

Wash water samples were - treated with ferric chloride at dosages of from 
150-190 ppm. Promising samples (150 and 170 ppm ferric chloride) were anal- 
yzed for fatty acid and total dissolved solids. While fatty acid analysis 
indicated that 90-97 percent of the acid was removed by the coagulation, 
total dissolved solids analysis showed that only 49-54 percent of the soap 
had been removed. The conclusion is that while coagulation is effective in 
removing fatty acid, it has little or no effect on the triethanolamine, glycer- 
ine, and perfume. No further pretreatment work has' been done with Neutrogena. 
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ML-ll BASED WASH WATER 


ML- 11 is another potassium fatty acid soap from Rochester Germicide 
Company. Our solids analysis indicated approximately 39 percent soap. Wash 
water prepared from this soap was treated with ferric chloride at dosages 
from 150-220 ppm (Table 10) . IR analysis indicates that ferric chloride at 
a dosage of between 200 and 220 ppm will remove 93-95 percent of the soap. 

Aluminum sulfate was also evaluated with ML-ll based wash water 
(Table 11) . As with other wash water systems investigated during the pro- 
gram, aluminum sulfate showed no promise as a primary coagulant. 

MIRANOL JEM BASED WASH WATER 

Miranol JEM based water was treated with various dosages of two coagu- 
lants: ferric chloride (Table 12) and aluminum sulfate (Table 13). In both 

instances the results were not promising. No further work was done with 
this cleansing agent. 

These sets of experiments conclude the chemical pretreatment portion of 
this effort. The overall results are discussed in Section V of this report. 
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TABLE 8 


Chemical Precipitation with Palmeto Soap 
Based Wash Water 


Aluminum Sulfate Series 


Experiment 

Number 

Coagulant 

FeCl 3 

(PPm) 

Floccula- 
tion Size 
and Rate 

(2) 

Settling 

Rate 

A-2132-A 

265 

1 

1 

A-2132-B 

250 

1 

1 

A-2132-C 

235 

1 

1 

A-2132-D 

220 

1 

1 

A-2132-E 

200 

1 

1 

A-2132-F 

350 

1 

1 

A-2132-G 

370 

1 

1 


Flocculation — rated on a scale of 0-5: 

0 - no floe 3 - faster growth 

1 - slow growth; small floe 4 - large floe in less than 1 min. 

2 - slow growth; larger floe 5 - large floe almost immediately 

than (1) 

Settling - rated on a scale of 0-5: 

0 - no floe 

1 - no settling 

2 - very little settling 

(less than 10%) 


3 

4 

5 


99-100% settled out in 20 min. 
99-100% settled out in 2-3 min. 
settles out in less than 2 min. 




TABLE 9 


Chemical Precipitation with Palmeto Soap Based Wash Water 
Aluminum Sulfate Series 


Experiment 

Number 

Aluminum 

Sulfate 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-3341-D-2 

160 

2 

1 (1) 

A-3341-C-2 

170 

2 

mm 

A-3341-B-2 

180 

2 

Hill 

A-3341 -A-2 

190 

2 

i < x > 


(1) Filtrate was cloudy 
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TABLE 10 


Chemical Precipitation with ML- 11 Soap 
Based Wash Water 

Ferric Chloride Series 


. 

Experiment 

Number 

\ 

Coagulant 

FeCl 3 

(ppm) 

Floccula- 
tion Size 
and Rate 

Settling 

Rate 

Residual 

Soap 

Content 

(ppm) 

Percent 

Soap 

Removed 

A-2136-A 

160 

0 

0 

- 

- 

A-2136-B 

170 

0 

0 

- 

- 

A-2136-C 

180 

5 

4 

- 

- 

A-2136-D 

190 

5 

4 

- 

- 

A-2136-E 

mm 

5 

4 

67 

93.3 

A-2136-G 

mm,/ 

5 

3 

66/63 

94.4/93.7 

A-2136-H 

210 

5 

3 

57/ - 

95.4/ - 

A-2136-I 

220 

5 

3 

48/53 

95.2/94.7 
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TABLE 11 


Chemical Precipitation with ML- 11 Soap Based Wash Water 
Aluminum Sulfate Series 


Experiment 

Number 

Aluminum 

Sulfate 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-2139-A 

160 

0 



A-2139-B 

180 

0 

0 

(1) 

A-2139-C 

200 

1 

0. 

(1) 

A-2139-D 

220 

2 

1 

(1) 

A-2139-E 

240 

2 

1 

(1) ' 

A-2139-F 

260 

3 

2 

(1) 


(1) Filtrate was cloudy 






























TABLE 12 ' 


Chemical Precipitation -with Miranol JEM Soap Based Wash Water 

Ferric Chloride Series 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-3340-A- 1 

120 

0 

0 

A-3340-B-1 

140 

0 

0 

A-3340-C-1 

160 

0 

0 

A-3340-D-1 

180 

0 

0 

A-3340-E-1 

200 

0 

0 

A-3340-F-1 

220 

0 

0 
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TABLE 13 


Chemical Precipitation with Miranol JEM Soap Based Wash Water 

Aluminum Sulfate Series 


Experiment 

Number 

Aluminum 

Sulfate 

.Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-3340-A-2 

120 

0 

0 

A-3340-B-2 

140 

0 

0 

A-3340-C-2 

160 

0 

0 

A-3340-D-2 

180 

0 

0 

A-3340-E-2 

200 

0 

0 

A-3340-E-2 

220 

0 

0 
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IV. DETERMINATION OF THE EFFECT OF 
ANTIFOAMS ON SOAP REMOVAL 


There was some concern at the beginning of the program that the use of 
antifoam agent with wash water would interfere with chemical precipitation of 
the soap. To determine the effect, if any, we mixed promising antifoams with 
synthetic wash water, ran chemical precipitation, and analyzed for residual 
soap. 

In our preliminary work under Contract NAS 9-14518, we found that 100 ppm 
of DB-110 had no adverse effect on the soap removal from an Ivory based wash 
water; and neither DB-110 nor SWS-211 - both at the 250 ppm level - had 

an adverse effect on removal with a Palmeto based wash water. 

During this program we determined the effect of other antifoams in com- 
bination with Ivory and Palmeto as well as ML-11. When 300 ppm of AF-75 was 
added to Palmeto based wash water there was no significant effect on the de- 
gree of soap removal; but as has been evident in other cases, a slightly lar- 
ger dosage of coagulant was required for effective soap removal (Table 14) . 

Antifoam AF-93 at a level of 200 ppm in ML-11 based wash water had a 
detrimental effect on the soap removal (Table 15) . Even with ferric chlor- 
ide dosages as high as 260 ppm, the soap removal was unsatisfactory. 

The presence of 200 ppm of SWS-211 had no adverse effect on soap removal 
from an ML-11 based wash water, however (Table 16) . Ferric chloride dosages 
of from 210 to 220 ppm were sufficient to remove 93-94 percent of the soap 
from the silicone-containing wash water. 

The removal performance was also satisfactory in the presence of 200 ppm 
of AF-75 antifoam; a dosage of 250 ppm ferric chloride was sufficient to re- 
move almost 96 percent of the soap from an ML-11 based wash water (Table 17) . 

The effect of antifoam on chemical precipitation was also determined 
for Ivory soap based wash water. In the presence of 200 ppm of DB-110 
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antifoam, a dosage of 180 ppm ferric chloride removed over 97 percent of the 
Ivory soap from the wash water (Table 18). 

The effect of AF-75 on removal was also negligible; a dosage of 180 ppm 
ferric chloride precipitated almost 94 percent of the Ivory soap from the 
wash water (Table 19) . 

Both antifoams AF-93 and AF-72 had a detrimental effect on soap removal 
when used in an Ivory based wash water. Ferric chloride dosages as high as 
260 ppm were ineffective in removing the soap (Tables 20 and 21) . 

In general, most of the silicone antifoams are compatible with the soap 
removal procedure, the two exceptions being; AF— 93 when used with Ivory or 
ML— 11, and AF-72 when used with Ivory. 
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TABLE 14 


Effect of AF-75^) Antifoam on Chemical Precipitation of 
Palmeto Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of 
Filtrate 
(ppm) 

Percent 

Soap 

Removed 

A-3341-A 

180 

5 

4 

- 

- 

A-3341-B 

190 

5 

4 

- 

- 

A-3341-C 


5 

4 

69/73 

93.1/92.7 

A-3341-D 

210 

5 

4 

• 83/86 

91.7/91.4 


(1) 300 ppm active solids 
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TABLE 15 


Effect of AF-93^) Antifoam on Chemical Precipitation of 
ML- 11 Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-3342-E-1 

170 

1 

2 (2) 

A-3342-F-1 

180 

2 

2 ( 1 2 ) 

A-3342-G-1 

190 

3-4 

2 (2) 

A- 3 342 -A- 1 

200 

4 

2 (2) 

A-3342-B-1 

210 

4 

2 ( 2 ) 

A-3342 -C-l 

220 

4 

2 ( 2 ) 

A-3342-D-1 

230 

4 

2 < 2 > 


240 

1 

2 

A-3347-B-1 

250 

1 

2 

A-3347-C-1 

260 

1 

2 


(1) 200 ppm active solids 

(2) Filtrate was cloudy 
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. TABLE 16 

Effect of SWS-21 1 U ) Antifoam on Chemical Precipitation 
of ML- 11 Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of 
Filtrate 
(ppm) 

Percent 

Soap 

Removed 

A- 3 343 -A- 1 

190 

4 

3 

- 

- 

A-3343-B-1 

200 

4 

3 

- 

- 

A-3343-C-1 

210 

4 

3 

61 

93.9 

A-3343-D-1 

220 

4 

4 

65/69 



(1) 200 ppm active solids 
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TABLE 17 


Effect of AF-75^) Antifoam on Chemical Precipitation of 
ML- 11 Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of 
Filtrate 
(ppm) 

Percent 

Soap 

Removed 

A-3343-E 

170 

0 

0 


- 

A-3343-F 

180 

1 

1 

- 

- 

A-3343-G 

190 

5 

2 

- 

- 

A-3342-A-2 

200 

5 

2 

- 

- 

A-3342-B-2 

210 

5 

2 

- 

- 

A-3342-C-2 

220 

5 

3 

- 

- 


230 

5 

. 3 

- 

- 

A-3347-A-2 

240 

5 

4 

- 

- 

A-3347-B-2 

250 

5 

4 

41 

95.9 

A-3347-C-2 

260 

5 

4 

- 



(1) 200 ppm. active solids 














































TABLE 18 


Effect of DB-llOU) Antifoam on Chemical Precipitation of 
Ivory Soap Based "Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of 
Filtrate 
(ppm) 

Percent 

Soap 

Removed 

A- 3 346 -A- 1 

180 

5 

3 

25 

97.5 

A-3346-B-1 

200 

5 

3 

- 

- 

A-3346-C- 1 

220 

5 

3 

- 

- 

A-3346-D-1 

240 

5 

3 

- 

- 


(1) 200 ppm active solids 
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TABLE 19 


Effect of AF-75^) Antifoam on Chemical Precipitation of 
Ivory Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of 
Filtrate 
(ppm) 

Percent 

Soap 

Removed 

A-3343-A-2 

180 

4 

3 

62 

93.8 

A-3343-B-2 ' 

200 

4 

3 

- 

- 

A-3343-C-2 

220 

4 

3 


- 

A-3343-D-2 

240 

4 

3 

- 

- 


(1) 20 0 ppm active solids 


























TABLE 20 


Effect of AF-72^) Antifoam on Chemical Precipitation of 
Ivory Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-3344-A-1 

170 

4 

3' ( 1 2 ) 

A-3344-B-1 

180 

4 

3 (2) 

A-3344-C-1 

190 

4 

3 (2) 

A-3344-D-1 

200 

4 

3 (2) 

A- 3 345 - A- 1 

220 

4 

2 (2) 

A-3345-B-1 

240 

4 

2 < 2 > 

A-3345-C-1 

- 260 

4 

2 (2) 


(1) 200 ppm active solids 

(2) Filtrate was cloudy 
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TABLE 21 


Effect of AF-93 ^ ) Antifoam on Cliemical Precipitation of 
Ivory Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccula- 

tion 

Size and 
Rate 

Settling 

Rate 

A-3344-A-2 

170 

3 

2 

(2) 

A-3344-B-2 

180 

3 

2 

(2) 

A-3344-C-2 

190 

3 

2 

(2) 

A-3344-D-2 

200 

3 

2 

(2) 

A-3345-D-1 

220 

2 

2 

(2) 

A-3345-E-1 

240 

2 

2 

(2) 

A-3345-F-1 

260 

2 

2 

(2) 


(1) 200 ppm active solids 

(2) Cloudy filtrate 
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V. SOAP AND ANTIFOAM SELECTION 


SOAP SELECTION 

A summary of the performance of wash waters based on various cleansing 
agents used throughout our study has been presented .in Table 22 , which indi- 
cates the amount of ferric chloride required for at least 90 percent soap re- 
moval, the quantity of antifoam necessary to defoam the given wash water, and 
the effect of the antifoam, if any, on the soap removal. The precipitation 
technique appears to work satisfactorily only with fatty acid soaps such as 
Ivory, Palmeto, and ML-11. The synthetics such as Neutrogena and Miranol JEM 
are not removed by the coagulants. 

As a group, the silicone antifoams are the most effective and are essen- 
tially identical in their performance with any given cleansing agent. With 
the exception of AF-72 and AF-93, the antifoams do not interfere with the 
soap removal. 

In order to keep the options open as to the physical form of the soap to 
be used, we recommended both a liquid and a bar soap. Unfortunately, only 
one bar soap gave satisfactory removal performance : that was Ivory . This 

soap can be defoamed with either AF-75 or DB-110, neither of which has a 
detrimental effect on the soap removal. 

There are two liquid soaps which gave good removal performance, but 
ML-11 has certain advantages over Palmeto. ML-11 is supplied as a 38 percent 
solution, twice as concentrated as the Palmeto. The use of ML-11 for space 
flight would minimize space and weight requirements for soap storage. 

Although ML-11 requires slightly larger dosages of ferric chloride for 
precipitation in the presence of antifoam (210 to 250 ppm versus 190 to 210 
ppm) , the antifoam dosage required is only one third of that, required for 
Palmeto . 
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'antifoam selection 

The anti foams to be used will he DB-110 with. Ivory soap and S.WS-211 
with ML— 11. in both cases these materials gave. the 'best. defoaming at the 
lowest dosage while permitting the highest degree of soap removal. 

In the third and fourth quarters of the 'program, we conducted studies 
to define operating parameters, using both Ivory and ML-11 soaps. All 
studies were conducted in the presence of the antifoams of choice (see the 
following Section VI - Optimization of Subsystem Parameters) . 
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TABLE 22 


Matrix Chart 

Comparative Performance of Wash Waters Based on Various Personal Cleansing Agents 



Optimum Ferric 

Promising Antifoams 

Effect of Antifoam on 

Soap 

Chloride Dosage 




Soap Removal 


ppm 

{Removal 

Efficiency) 

Reference 

Number 

Material 

Dosage 

(ppm) 

Reference 

Number 

Effect 

FeCl 3 

Dosage 

(ppm) 

Reference 

Number 




AF-75 

300 

A- 2 140 

Negligible 

200-210 

A-3341 

Palmeto 

175-180 

A-3337 

SWS-211 

300 

A-2140 

Negligible 

200-210 

A-3338 


(92-95%) 


DB-110 

300 

A-2140 

Negligible 

190 

A-3337 




AF-93 

100-200 

A-2142 

Destroys 

Removal 


A- 3 344, -5 


170-185 

A-2114 

AF-75 

200 

A-2142 

Negligible 

180 

A-3343 

Ivory 

(95-97%) 

AF-72 

200 

A-2142 

Destroys 

Removal 

- 

A-3344i -5 




DB-110 

200 

A-2142 , 

Negligible 

180 

A-3346 

Neutrogena 

Ineffective 
(approx. 50%) 

A-2137 

NA< 1 ' 2 > 

NA 

A-2143 

NA (1 > 

NA 

NA 




AF-93 

100 

A-2141 

Destroys 


A-3342,47 


200-220 

A-2136 

Removal 



ML- 11 

(93-95%) 

AF-75 

200 

A-2141 

Negligible 

250 

A- 3 347 



SWS-211 

100 

A-2141 

Negligible 

210-220 

A-3343 




AF-75 

50 

A-2143 

NA 

- 


Miranol 

Ineffective 

A-3340 

DB-110 

50 

A-2143 

NA 

- 

- 

JEM 



Colloid 680 

200 

A-2143 

NA 

~ 

- 


(1) NA - Not Applicable 

(2) Neutrogena -based wash water alone gave negligible foam 




VI. OPTIMIZATION OF SUBSYSTEM PARAMETERS 


Under contract NAS 9-14518, and during the first two quarters of the 
present program we investigated soap removal agents and antifoams for vari- 
ous candidate personal cleansing agents. The culmination of this work was 
the selection of two soaps, one liquid and one bar soap, and two antifoam 
agents for further study. During the last two quarters our efforts were 
directed toward optimizing the subsystem parameters using both of the can- 
didate soaps and the antifoams. 

The two soaps selected were Ivory bar soap and ML-11 liquid soap. To 
wash waters prepared from these soaps, we added DB-110 and SWS-211 silicone 
antifoams , respectively . 

REPRODUCIBILITY 

One important design consideration is the reproducibility of the re- 
moval procedure. To investigate this we conducted identical coagulation/ 
flocculation jar experiments in triplicate using identical ferric chloride 
dosages for each jar; the evaluations were conducted on synthetic wash 
water containing 0.10 percent soap solids. The experiments were con- 
ducted using a gang stirrer and dasher mixer as previously described. 

The results for both Ivory Soap and ML-11 are summarized in Table 23. 

Within the limits of accuracy of the analysis procedure, the removal 
technique is reproducible for both soaps ; reproducibility should present 
no problem in the design and operation of a removal system based on this 
concept. 

EFFECT OF SOAP CONCENTRATION IN 
THE WASH WATER ON THE REMOVAL 

To determine the effect of concentration, we prepared synthetic wash 
waters from both candidate soaps at concentrations of 0.05 percent, 0-10 
percent, and 0.25 percent. Coagulation and flocculation studies were again 
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conducted using the gang stirrer. The ferric chloride dosages for each 
jar were adjusted in proportion to the soap concentration. At each soap 
concentration we investigated several ferric chloride dosages over a nar- 
row range in the vicinity of what was estimated to be the "optimum" dosage. 
Promising samples were filtered and analyzed for residual soap using the 
previously described carbon tetrachloride extraction/IR technique. 

The results of these investigations are presented in Tables 24 and 25 
for the two soaps. In both cases the residual soap content in the fil- 
trate of the promising samples increased with increasing concentration of 
soap in the synthetic wash water feed. The increase was not proportional, 
however, and as a result the percent soap removed tended to increase with 
increasing soap content-of the wash water feed.’ 

The variations in removal are not significant and in all cases the re- 
moval is in excess of 90 percent. This should present no problem in the 
design and operation of a pretreatment system based on this concept since 
the trace amount of residual soap will be removed by carbon adsorption. 

EFFECT OF THE DURATION OF SLOW STIRRING 
ON THE DEGREE OF SOAP REMOVAL 

In conducting our coagulation and flocculation studies over the course 
of the last two programs we have standardized on a fast mix time of 30 sec- 
onds, and a slow mix time of 15 minutes, a time commonly employed in jar 
water clarification studies. 

In order to determine the optimum slow mix time we ran series of iden- 
tical jar experiments. Soap concentration and ferric chloride dosage were 
kept constant. Samples were mixed for periods of 0, 2, 5, 15, 30, and 60 
minutes, filtered, and analyzed for residual soap content. The results 
are presented in Tables 26 and 27. For slow mix periods of 0 and 2 min- 
utes the experiments were done in duplicate to provide greater reliability. 

With Ivory soap based wash water, maximum soap removal was attained 
with no slow mixing; with ML-11 based wash water there appears to be a 
fall-off in soap removal with slow mix periods of less than 5 minutes. 
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There is not a good explanation as to why those samples mixed for 
2 minutes had poorer soap removal than those with no slow mixing. How- 
ever, this discrepancy will probably be resolved when the optimum slow mix 
period is determined as part of the Bread Board demonstration under con- 
tract NAS 9-15369. 

EFFECT OF HIGHER FERRIC CHLORIDE 
SOLUTION CONCENTRATION 

For convenience in all of our coagulation and flocculation work we 
have used 1 percent solutions of all coagulants. In order to conserve 
weight and space during missions of long duration, however, it will be 
necessary to employ more concentrated solutions. 

To investigate the effect of higher concentration of ferric chloride 
on soap removal, we employed a 10 percent solution, and dosages in the 
range which had already proved to be near optimum. The results are presen- 
ted in Tables 28 and 29. 

For Ivory soap based wash water, a dosage of 175 ppm of ferric chloride 
applied from a 10 percent solution removed over 98 percent of the soap. 

This result compares favorably' with previous removal results when using 
1 percent solution. 

For ML-11 based wash water, a dosage in the range of 230 ppm resulted 
in the removal of approximately 93.5 percent of the soap, comparable with 
previous results using 1 percent ferric chloride solution. 

EFFECT OF FILTER PORE SIZE AND OPERATING 
PRESSURE ON THE DEGREE OF SOAP REMOVAL, 

WATER FLUX, AND RATE OF FLUX DECLINE 

To evaluate filtration characteristics, we employed a cartridge type 
pressure filtration housing. Type 1-M {AMF Cuno Division) . This housing 
holds single cylindrical cartridges approximately 25 cm long by 6.7 cm in 
diameter, supplied in various pore sizes. 

In order to process large quantities of water, the filter housing was 
mounted to the bottom flange of a 38-liter glass-lined reactor, used here 
as a reservoir and pressure vessel (Figure 1) . 
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Wash water based on-ML-11 soap was prepared in 15 Kg batches in the 
laboratory and was coagulated with ferric chloride at a dosage of 240 ppm 
using a laboratory stirrer. Coagulated wash water was charged to the re- 
actor through the fill pipe at the top. Nitrogen pressure was applied to 
the reactor through the top, and water flow rate was controlled by the 
shut-off valve downstream from the filter housing. To keep the soap sludge 
suspended, the reactor was stirred for 30 seconds before each filtration 
experiment. 

To measure the flow rate through the filter, samples of water were col- 
lected at the desired pressure for timed intervals. The water was weighed 
and volume was computed from the weight of water. 

Since the final system will function on a batch basis and will operate 
intermittently, we followed the flow rate decline for the filters as a func- 
tion of cumulative volume rather than time. The tests were terminated 
when the flow rate fell below approximately 0.5 liters per minute. Water 
aliquots were sampled occasionally for analysis to determine the residual 
soap content. Since the filter housing was constructed from clear plastic 
we were also able to observe the filtering action and note any channeling 
or plugging in the housing. 

The filtrate was also observed for cloudiness or obvious signs of 
sludge. 

Early microscopy work on soapy wash water coagulated with ferric chlor- 
ide revealed large soap aggregates, 500 to 750 microns in diameter, with 
smaller fragments measuring 50 to 75 microns. Therefore, our first filtra- 
tion experiment employed a nominal 50-micron pore depth filter {Micro-Klean 
G78L8 - AMF, Cuno Division) . The flow rate versus cumulative volume data 
for this filter is presented in Table 30. 

The 50-micron filter allowed a large flow rate - nearly 28 liters per 
minute at 5 psig - but failed to remove most of the sludge from the efflu- 
ent, even after processing over 40 liters of water. 

A closer examination of the sludge under a microscope revealed that the 
relatively large soap aggregates were actually composed of 5 to 20 micron 
particles. The aggregates appear to break apart under the relatively high 
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shear associated with the filtering operation, with the small particles 
passing through coarse filters. 

Next, a nominal 10-micron pore depth filter was evaluated (Micro- 
Klean G78C8 - AMP, Cuno Division). The initial flow rate was high, with a 
fairly clear effluent for the first 10 to 12 liters. The residual soap con- 
’ tent of the first effluent was high - 94.5 percent - but had dropped' to 
90.5 percent by the time 13 liters of wash water had been processed (Table 
31) . By the time 20 liters of water had been filtered, large quantities of 
sludge were coming through the filter. Because of the poor filtering, the 
test was terminated after a cumulative volume of 30 liters. 

A 5-micron depth filter (Micro-Klean G78B3 - AMF, Cuno Division) was 
evaluated next at two pressures: 5 and 25 psig. At 5 psig the initial 

flow rate was again high, with over 94 percent of the soap being removed 
from the effluent (Table 32) . After 16 liters of water had been treated, 
however, the filtrate became slightly cloudy? soap analysis indicated that 
the removal had dropped to less than 88 percent. 

With subsequent filtering, the cartridge began to blind and the efflu- 
ent cleared. By the time 25 liters of water had been recovered, the soap 
removal had returned to over 94 percent. 

At 25 psig the 5-micron filter permitted a tremendous flow rate, in 
excess of 60 liters per minute. Large amounts of sludge came through the 
filter from the beginning of the test, which was terminated after 22 liters 
of water- had been processed. 

A 1-micron depth filter (Micro-Wynd DCCPY - AMP, Cuno Division) was 
evaluated at two pressures - 10 and 25 psig - and gave consistently high 
soap removal throughout the evaluations. 

A key criteria in evaluating filters is the decline in flow rate ver- 
sus cumulative volume processed. This information has been presented 
graphically in Figure 2, as plots of flow rate in liters per minute versus 
cumulative liters processed. The curves were constructed from the data 
in Tables 31, 32, 33, and 34. 


- 40 - 



As was expected, the filters with larger pore diameters allowed 
higher flow rates, both initially and at low cumulative volumes (i.e., 

10 to 30 liters) . Because of the relatively large quantity of suspended 
matter in the water, however, all the filters, regardless of pore size, 
blinded with suspended matter. As a result, the flow rates tended to con- 
verge and approach asymptotically some minimal flow rate. This rate cor- 
responds to that for a "sludge cake" filter, which of course could be the 
same in each case, regardless of filter type. 

Despite the rapid drop in flow rate for all of these filters, the 
rate remained at an acceptable level - approximately 500 cc per minute for 
a cumulative volume of greater than 60 liters. This volume of water is 
equivalent to the requirements for a four-man crew for hand washing for a 
week assuming a nominal 40 washes per day and 225 cc per wash (0.5 pound 
per wash) . 

It should be noted that at no point during the testing did the sludge 
collected by the filters exceed the capacity of the filter housing, al- 
though in the case of the 1-micron filter evaluated at 25 psig with 64 
liters of water, the sludge had filled to the wall of the housing in the 
lower 1/4 of the housing. 

A cursory evaluation indicated that at 5 psig the solids content of 
the sludge was approximately 12 percent; this is equivalent to a potential 
water recovery in excess of 99 percent for a pretreatment system based on 
coagulation and flocculation followed by pressure filtration. 

To summarize the results of the filtration tests, all of the filters 
with the exception of the 1-micron filter are too large in pore .size and 
permit at least some sludge to come through into the effluent. The drop 
in flow rate with cumulative volume is rather severe for all filters inves- 
tigated. Increased pressure does not have a dramatic effect on the abso- 
lute flow rate but does permit a greater cumulative volume before an unac- 
ceptably low rate is reached {e.g. , 500 cc/minute) . A comparison of 
Tables 33 and 34 illustrates the difference. With 25 psig pressure, 50.3 
liters of water was processed before a flow rate of 0.43 liters per minute * 
was achieved, while at 10 psig only 30.1 liters was filtered before a flow 
rate of 0.42 liters per minute was reached. 
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TABLE 23 


Determination of -the Reproducibility of the Soap Removal 
Using Ferric Chloride Coagulant 


Experiment 

Number 

Flocculation 
Size and Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 

Ivory Soap Based Wash Water 

with ; 200 

ppm DB 110 Antifoam, 

0.10 Percent Soap; 180 ppm Ferric Chloride Dosage 


A- 3349 -A 

5 

4 

25/25 

97.5/97.5 

A-3349-B 

5 

4 

25/25 

97.5/97.5 

A-3349-C 

5 

4 

24 

97.6 

ML-11 Soap Based Wash Water 

with: 0.10 Percent Soap, 


200 ppm SWS- 

-211 Antifoam; 220 ppm Ferric Chloride Dosage 

A-3347-C 

5 

4 

45/49 

95.5/95.1 

A-3347-B 

5 

4 

57 

94.3 

A-3347-A 

5 

4 

56/57 

94.4/94.3 
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TABLE 24 


The Effect of Soap Concentration in the Synthetic Wash Water 
on the Degree of Soap Removal - Ivory Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccu- 
lation 
Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 

0.05 Percent Soap in the Wash Water 



A-2148-D 

85 

4 

4 

15/17 

97.0/96.6 

A-2148-E 

90 

4 

4 

16/12 

96.8/97.6 

0.10 Percent Soap in the Wash Water 



A-2114-L 

175 

5 

5 

37 

96.3 

A-2114-G 

185 

5 

5 

29 

97.1 

0.25 Percent Soap in the Wash Water 



A— 2148-A 

425 

5 


45 

98.2 

A-2148-B 

440 

5 

- 

39/41 

98.4/98.4 
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TABLE 25 


The Effect of Soap Concentration in the Synthetic Wash Water 
on the Degree of Soap Removal - ML-11 Soap Based Wash Water 


Experiment 

Number 

Ferric 

Chloride 

Dosage 

(ppm) 

Floccu- 
lation 
Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 

0.05 Percent Soap in 

the Wash Water 



A— 2147-A 

100 

4 

1 

29/35 

94.2/93.0 

A-2147-B 

105 

4 

1 

36 

92.8 

0.10 Percent Soap in 

the Wash Water 



A-2136-H 

210 

5 

3 

57 

94.3 

A-2136-I 

220 

5 

3 

48/53 

95.2/94.7 

0.25 Percent Soap in 

the Wash Water 



A-2146-A 

500 

5 

. 3 

85/91 

96.6/96.4 

A-2146-B 

525 

5 

3 

72/71 ' 

97.1/97.2 














TABLE 26 


The Effect of Duration of Slow Stirring on the Degree of Soap Removal 

Ivory Soap Based Wash Water - 0.10 Percent Soap, With 
200 ppm DB-110 Antifoam, 180 ppm Ferric Chloride 


Experiment 

Number 

Duration - of 
Slow Mixing 
{minutes) 

Floccu- 
lation 
Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 


0 

4 

4 

18/18 

98*2/98.2 

A-6053-B 

0 

4 

4 

19/19 

98.1/98.1 

A-6053-C 

2 

4 

4 

19/21 

97.9/98.1 

A— 605 3-D 

2 

4 

4 

19/19 

98.1/98.1 

A-3349-D 

5 

5 

4 

25/25 

97.5/97.5 

A— 3349— E 

15 

5 

4 

24 

97.6 

A-3349-F 

30 

5 

4 

21/22 

97.9/97.8 

A-3349-G 

60 

5 

4 

25 

97.5 
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TABLE 27 


The Effect of Duration of Slow Stirring on the Degree of Soap Removal 

ML-11 Soap Based Wash Water - 0.10 Percent Soap With 
200 ppm SWS-211 Antifoam, 220 ppm Ferric Chloride Dosage 


Experiment 

Number 

Duration of 
Slow Mixing 
(minutes) 

Floccula- 
tion 
Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 

A-6062-A 

0 

5 

- 

62 

93.8 

A-6062-B 

0 

5 

- 

62/62 

93.8/93.7 

A-6062-C 

2 

5 

- 

91/83 

90.9/91.7 

A-6062-D 

2 

5 

- 

84/83 

91.6/91.7 

A-3348-A 

5 

5 

4 

45/55 

95.5/94.5 

A— 3348-B 

15 

5 

4 

47/50 

95.3/95.0 

A-3348-C 

30 

5 

4 

49/48 

95.1/95.2 

A— 3348-D 

60 

5 

4 

38/37 

96.2/96.3 
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■TABLE 28 


The Effect of Ferric Chloride Solution Concentration 
on the Degree of Soap Removal 

Ivory Soap Based Wash Water with 200 ppm DB-110 Antifoam 


Experiment 

Number 

Ferric Chloride 
Dosage from 10% 
Solution (ppm) 

Floccula- 
tion 
Size and 
Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 

A-3350-A 

175 

5 

3 

16/19 

98.4/98.1 

A-3350-B ' 

180 

5 

3 

27/28 

97.3/97.2 

A-3350-C 

185 

5 

3 

- 

- 

A-3350-D 

190 

5 

3 

- 

- 
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TABLE 29 


‘The Effect of Ferric Chloride Solution Concentration 
on the Degree of Soap Removal 

ML— 11 Soap Based Wash Water with 200 ppm SWS-211 antifoam 


Experiment 

Number 

Ferric Chloride 
Dosage from 10% 
Solution (ppm) 

Floccula- 

tion 

Size and 
.Rate 

Settling 

Rate 

Actual Soap 
Content of the 
Filtrate (ppm) 

Percent Soap 
Removed 

A-6062-E 

210 

5 

' ' 4 


- 

A-6062-F 

220 

5 

4 

- 

- 

A-6062-G 

230 

5 

4 

64/66 

93.6/93.4 

A-6062-H 

240 

5 

4 

57 

94.3 
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TABLE 30 


Flow Rate Vs. Cumulative Volume for Pressure Filtration of 
Coagulated Wash Water Through a 50-Micron Filter at 5 psig ^ 


Cumulative Volume 
of Water, Liters 
(Gallons) 

Flow Rate, 
Liters/Min. 
(Gal. /Min. ) 

Residual 
Soap Content 
(ppm) 

Percent 

Soap 

Removed 

9.2 (2.4) 

27.6 (7.29) 

(2) 

- _ 

17.7 (4.7) 

25.4 (6.7) 

(2) 

- 

28.8 (7.6) 

22.3 (5.9) 

(2) 

- 

40.5 (10.7) 

23.4 (6.2) 

(2) 

- 


(1) 0.1 Percent ML-11 soap based wash water 
coagulated with ferric chloride. 

(2) Large quantities of soap/ferric chloride 
sludge came through the filter. 
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TABLE 31 


Flow Rate Vs. Cumulative Volume for Pressure Filtration of 
Coagulated Wash Water Through a 10-Micron Filter at 5 psig 


Cumulative Volume 
of Water, Liters 
(Gallons) 

Flow Rate, 
Liters/Min. 
(Gal. /Min. ) 

Residual 
Soap Content 
(ppm) 

Percent 

Soap 

Removed 

mmm 

28.8 (7.6) 

55 

94.5 

13.3 (3.5) 

24.5' (6.5) 

93/97 

90 .7/90. 3 

20.9 (5.5) 

15.1 (4.0) 

(3) 

- 

25.8 (6.8) 

9.9 (2.6) 

(3) 

- 

30.3 (8.0) 

9.0 (2.4) 

(3) 

- 


(1) 0.1 Percent ML-11 soap based wash water coagulated 
with ferric chloride. 

(2) Small quantities of soap/ferric chloride sludge 
came through the filter. 

(3) Large quantities of soap/ferric chloride sludge 
came through the filter. 
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TABLE 32 


Flow Rate Vs. Cumulative Volume for Pressure Filtration of 
Coagulated Wash Water Through a 5 -Micron Filter at 5 psig 


Cumulative Volume 
of Water, Liters 
(Gallons) 

6.8 

(1.8) 

11.5 

(3.0) 

16.3 

(4.3) 

20.9 

(5.5) 

24.9 

(6.6) 

31.0 

(8.2) 

36.9 

(9.7) 

41.9 

(11.1) 

48.9 

(12.9) 

55.1 

(14.6) 

57.7 

(15.2) 


Flow Rate, 
Liters/Min . 
(Gal. /Min. ) 


27.2 (7.2) 


18.7 (4.9) 


2.6 


4.5 (1.2) 


4.0 (1.1) 


.82 


(0.40) 


.26 


0.7 (0.18) 


1.0 (0.26) 


0.5 (0.13) 


Residual I Percent 

Soap Content 

(ppm) I Removed 


57/59 


120/127 ' ' | 88.0/87.3 


55/56 


94.5/94.4 


0.1 Percent ML-11 soap based wash water coagulated 
with ferric chloride. 


Filtrate was slightly cloudy. 
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TABLE 33 


Flow Rate Vs. Cumulative Volume for Pressure Filtration, of 
. Coagulated Wash Water Through a 1-Micron Filter at 10 psig 


( 1 ) 


Cumulative Volume 
of Water, Liters 
(Gallons) 

Flow Rate, 
Liters/Min. 
(Gal. /Min. ) 

Residual 
Soap Content 
(ppm) 

Percent 

Soap 

Removed 

4-7 

(1-2) 

19.0 

(5.0) 

42/42 

95.8/95.8 

mm 

(1.9) 

10.0 

(2.6) 

- 

- 

9.2 

(2.4) 

3.8 

(1.0) 

- 

- 

12.0 

(3.2) 

2.8 

(0.74) 

- 

- 

16.3 

(4.3) 



39/39 

96.1/96.1 

21.0 

(5.5) 


(0.42) 

- 

- 

24.5 

(6.5) 

IB 


- 

- 

27.6 

(7.3) 

0.61 

(0.16) 

- 

- 

30.1 

(8.0) 

0.42 

(0.11) 

- 

- 

33.8 

(8.9) 

0.41 

(0.11) 

43/47 

95.6/95.3 


(1) 0.1 Percent ML-11 soap based wash water coagulated 

with ferric chloride. 


- 52 - 
































TABLE 34 


Flow Rate Vs. Cumulative Volume for Pressure Filtration of 
Coagulated Wash Water Through a 1-Micron Filter at 25 psig 


Cumulative Volume 
of Water, Liters 
(Gallons) 

Flow Rate, 
Liters/Min. 
(Gal. /Min. ) 

Residual 
Soap Content 
(ppm) 

Percent 

Soap 

Removed 

7.2 

(1.9) 

28.8 

(7.6) 

43 

95.7 

10.1 

(2.7) 

EB 


- 

- 

13.3 

(3.5) 

6.3 

(1.7) 

- 

- 

16.9 

(4.5) 

3.6 

(0.95) 

- 

- 

23.0 

(6.1) 

2.0 

(0.53) 

- 

- 

28.4 

(7.5) 

n 


- 

- 

33.5 

(8.8) 

0.73 

(0.19) 

- 

- 

38.1 

(10.0) 

0.57 

(0.15) 

- 

- 

44.7 

(11.8) 

0.47 

(0.12) 

- 

- 

50.3 

(13.3) 

0.43 

(0.11) 

- 

- 

57.5 

(15.2) 



- 

- 

64.2 

(17.0) 

0.33 

(0.09) 

44 

95.6 


(1) 0.1 Percent ML-11 soap based wash water coagulated 

with ferric chloride. 



















' FIGURE 1 

Pressure Filtration Apparatus 
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FIGURE 2 


Flow Rate Vs. Cumulative Volume of Water Processed 
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TRADE-OFF ANALYSIS 


While this trade-off analysis in a strict sense applies to the present 
contract, it is really the culmination of three programs covering a four- 
year period leading up to the design and construction of a breadboard wash 
water renovation system. Where our analysis bears on previous work the ap- 
propriate contracts will be referenced. 

The bulk of this analysis will be substantiated by the body of the final 
report, which has preceded this section. Again, references to pertinent 
sections will be made to support the analysis. 

SOAPS, ANTIFOAMS, AND COAGULANTS 

Numerous different cleansing agents have been evaluated as part of the 
three development programs. These materials can be categorized as either 
fatty acid soaps or synthetics, as follows: 

. Synthetics : 

Miranol JEM (The Miranol Company) 

Neutrogena (Neutrogena Corporation) 

Igepon CN-42 (GAF Corporation) 

Aerosol 18 (American Cyanamid) 

Aerosol MA-80 (American Cyanamid) 

. Fatty Acid Soaps : 

Olive Leaf Soap - discontinued (Rochester Germicide) 

Palmeto (Rochester Germicide) 

ML- 11 (Rochester Germicide) 

Ivory Soap 

The cleansing agents were selected for evaluation on the basis of their 
apparently good dermatological characteristics; the three main criteria for 
final selection of a cleansing agent were ability to.be removed from wash 
water by chemical precipitation, ability to be defoamed using antifoam 
agents, and compatibility of the required antifoam agent with the chemical 
precipitation process. 
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Under contracts NAS 9-14518 and NAS 9-14965 a study was conducted to 
identify coagulants for the candidate soaps. These studies revealed that 
only the fatty acid soaps are successfully removed by chemical precipitation, 
while the synthetics are generally unaffected. Aluminum sulfate was sur- 
prisingly unsuccessful with the various cleansing agents, but ferric chlor- 
ide has proved to be effective with all of the candidate fatty acid soaps. 
This agent will be used for all future work. 


As part of the coagulation work, several polymeric coagulants includ- 
ing anionics, cationics, and nonionics were evaluated in conjunction with 
primary coagulants. None of these materials enhanced the soap removal and 
therefore will not be employed in future work. 


Under contracts NAS 9-14518 and NAS 9-14965 antifoam agents were evalu- 
ated for wash waters based on each of the candidate soaps. Antifoams inves- 
tigated can be categorized as follows: 


Material 


Manufacturer 


Type 


Nalco 73C44, 73C37 

Nalco 71-D5, 8616, 867, 
and 4WP126 

Dow Corning DB-110 

Antifoam 71, AF-72, AF- 
75, AF-93 

SWS-211 

Foamaster, Foamaster B, 
VL, AP, G, and TMC 

Foamex/Pegosperse 10 OL 

Colloid 680, 68 IF 

Nalco 121 


Nalco Chemical Company 
Nalco Chemical company 
Dow Corning Corporation 


Amides 

Fatty acid 
esters 

Silicone 


General Electric Company Silicone 


SWS Silicones Corporation 

Nopco Chemical Division, 
Diamond-Shamrock 

Glyco Chemical 

Colloids, Inc. 

Nalco Chemical Company 


Silicone 

Proprietaries 

Proprietary 

Proprietary 

Proprietary 


Ferric chloride, cal- 
cium chloride, aluminum 
sulfate, and magnesium 
sulfate 


Heavy metal 
ions 
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Based on both, a NASA-JSC "shake test" and a dynamic test using cas- 
cadingwater, the silicones as a group were far superior to any of the 
other types of antifoam agents. 

The rationale for final soap and antifoam agent selection can be 
found on pages 33 and 34 of this report, and in Table 22 on page 35, which 
is a matrix chart comparing the soap removal performance of wash waters 
based on various personal cleansing agents when used in combination with 
the more promising antifoams. 


NASA-JSC requested that the selection include one liquid and one bar 
soap; the choices of soap and respective antifoam agent as well as required 
dosage of ferric chloride for optimum soap removal are as follows: 


Coagulant Dosage 
Required 

Antifoam Agent (ppm Ferric 
Soap Type Soap (dosage, ppm) Chloride) 


Liquid ML-11 

Bar soap Ivory 


SWS-211 (100) 210-220 

DB-110 (200) 180 


(1) Assumes a wash water based on 0.1 percent soap 


SOAP CONCENTRATION 

In a final pretreatment system, it is possible that there would be 
relatively small variations in the soap concentration of the wash water to 
be treated. 

Within reasonable limits, the percentage of soap removed from wash 
water is unaffected by the concentration of soap in the wash water when 
chemical precipitation is used - assuming that the dosage of ferric chlor- 
ide is in proper proportion to the soap. The removal performance for soap 
concentrations in the range of from 0.05 percent to 0.25 percent, for both 
Ivory and ML-11 soaps, is presented on page 36 of this report, and in 
tables 24 and 25 on pages 43 and 44, respectively. 
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COAGULANT CONCENTRATION 


In a final system it may be necessary to employ the coagulant in a 
more concentrated form than has been used -in the laboratory in order to 
conserve weight and space. Depending on the quantities of water to be 
treated, the concentration may have to be as high as 10 percent. 

During the current program wash water was treated with 10 percent 
ferric chloride. There was no adverse effect on removal, with results 
being essentially the same as those obtained with 1 percent ferric chlor- 
ide solutions. For details see page 38 of this report, and tables 28 and 
29, pages 47 and 48, respectively. 

DURATION OF MIXING 

In a final system, to minimize the cycle time and therefore reduce 
turn-around time for the equipment as well as power requirements, it would 
be desirable to have as short a mix time as possible after the addition of 
the ferric chloride coagulant. 

An investigation of mixing time spanning from "no mixing" up to 60 min- 
utes revealed that mixing for periods longer than 5 minutes does not en- 
hance the soap removal. Mixing time of less than 5 minutes (i.e., 0 to 
2 minutes) had no adverse effect on the soap removal with Ivory soap based 
wash water; but with ML-11 soap based wash water, the removal was somewhat 
reduced. The data as presented in tables 26 and 27 is not completely con- 
vincing for ML-11 based wash water, and a .repeat of this investigation 
using a breadboard model will be required to further- define the need for 
slow mixing. 

SOAP REMOVAL 

In order to remove the coagulated soap from the wash water, it will 
be necessary 'to use some form of filtration. There are numerous types of 
filtration systems for industrial use - including rotary vacuum, plate- 
and- frame pressure filtration, batch leaf, regenerable or backf lushing , or 
in-line cartridge filtration using either pleated membrane or depth filter 
type cylindrical cartridges. Many of these systems - such as rotary vacuum. 


- 59 - 



plate-and-f rame , and batch leaf filters - have large space weight and energy 
requirements and are overdesigned for the present application. 

Regenerable filters are . attractive because of their low maintenance 
requirements, but because of the relatively high level of suspended matter 
in the coagulated wash water, large quantities of water would have to be 
used for frequent backf lushing. 

A simplified pressure filter or in-line cartridge filter is the most' 
promising for the present application. Coagulation of an Olive Leaf Soap 
based wash water with 170 ppm ferric chloride followed by pressure filtra- 
tion using an asbestos pad depth filter removed in excess of 95 percent of 
the soap from the water in investigations conducted under contract NAS 9- 
13536. 

THE EFFECT OF POKE SIZE ON REMOVAL 

Investigations were conducted on depth type filter cartridges employing 
media with nominal pores in the range of from 50 down to 1 micron. Soap 
analysis and visual results indicate that only the 1-micron filter provides 
for a consistently high level of soap sludge removal when using both high 

i 

and low pressures. -Filtration work has been discussed on pages 38 through 
41. 

WATER FLOW RATE AND FLOW RATE DECLINE 

The initial flow rates for the various pore sized filter cartridges de- 
pend, as might be expected, on the nominal size of the pores. During filtra 
tion the rate drops off rapidly, however, as a result of blinding of the 
filters caused by the relatively high level of suspended matter in the coagu 
la ted wash water. 

The 1-micron filters are the only ones capable of consistently removing 
all of the sludge and therefore are the only cartridges of interest here. 

Once a cumulative volume of 30-40 liters of water has been processed 
with these filters, there is little difference in the flow rate between the 
cartridge when operated at 10 psig versus a 25 psig pressure. Increased 
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pressure, however, does offer the advantage of allowing a greater cumula- 
tive volume of water to be filtered before the flow rate reaches an unac- 
ceptably low level. 

For a final system, a staging of pressures might be desirable, begin- 
ning with a 10-15 psig pressure and gradually increasing the pressure to 
50 or 100 psig as the sludge content of the filter housing reaches near 
capacity. This would permit higher cumulative amounts of water to be pro- 
cessed before a low flow rate required a replacement of the unit. 

Based on the work conducted under the present program (see pages 38 
through 41, we feel that a 1-micron cartridge, when used in a pretreatment 
system, would be capable of processing a minimum of one week’s requirement 
of hand wash water for a four-man crew. Rather than rely on the replacement 
of filter cartridges for such a pretreatment system, we envision the use of 
several housings in parallel, each containing a fresh cartridge. As one 
housing became filled, another would be valved into the system to take its 
place. 

Housings are available which can accommodate two or three cylinderical 
cartridges stacked end to end. Based on our present work, we would esti- 
mate that a cluster of eight housings, each containing three cartridges, 
would be able to process all of the pretreated hand wash water generated 
by a four-man crew during a six-month mission. 

THE CONCEPT FOR A PRETREATMENT SYSTEM 

We envision that a typical pretreatment system based on the chemical 
precipitation concept would be as follows : 

Water from a hand washing station would be picked up and passed through 
a prefilter to remove hair, lint, and other suspended matter. The water 
sould comprise either Ivory or ML-11 fatty acid soap in a concentration of 
•from 0.05 - 0.25 percent as well as dissolved salts at approximately 150 
ppm consisting primarily of the sodium salts of chloride and sulfate. 

Following the prefilter, antifoam would be added to the water at the 
appropriate dosage (e.g., 200 .ppm for a 0.1 percent solution of Ivory soap) 
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using a metering pump. The antifoams , DB-110 with Ivory soap, and SWS-211 
with ML-11 , would be dispensed from a nominal 4 percent dispersion in water. 
This concentration provides for a reasonably stable dispersion and control 
of the dosage amount in the range of interest. 

Once the antifoam has been added, the water would enter a centrifugal 
liquid/air separator. This separator would likely be a small rotating vane 
type separator such as the one developed by Hamilton-Standard for the Sky 
Lab program for urine separation. The antifoam in the water should aid in 
preventing excessive foaming in the separator. 

Water exiting from the separator would be picked up by centrifugal 
pump, which would carry it to a mixing chamber. Prior to the mixing cham- 
ber, ferric chloride, the coagulating agent, would be added to the water 
from a nominal 10 percent solution - again using a metering pump . The do- 
sage would depend on the soap being used and the concentration of that soap 
in the wash water, but would typically be in the range of from 100-600 ppm. 
The water, now dosed with ferric chloride, would enter a mixing chamber, 
where it would be slow-mixed at 25-50 rpm for a period of 2-5 minutes. 

This mixing chamber would consist of a nominal 10 cm diameter cylinder with 
a free-floating piston inside to allow an adjustable volume inside the cav- 
ity. At the bottom of the chamber would be an impeller to gently agitate 
the water. 

Once the water had slow-mixed, the piston would be driven downward by 
gas pressure to drive the water out of the cylinder and through a depth 
type filter cartridge. This cartridge would be a nominal 1-micron pore 
size filter such as the Micro-Wynd DCCPY (AMP, Cuno Division) or. suitable 
substitute, inside a conventional cartridge housing such as AMF-Cuno Unit 
Type IB. 

This would comprise the entire pretreatment system; water from' this 
unit might be further processed by multifiltration through a carbon adsor- 
ber and deionizer, through a reverse-osmosis unit, through an electrodi- 
alysis apparatus or some other combination of the above to remove traces 
of organics and electrolytes. 
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